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Vibration of an Elastic Globe with a Homogeneous 
Mantle over a Homogeneous Core. 
Vibrations of the First Class 


By 
Yasuo SATO and Tosimatu MATUMOTO, 


Earthquake Research Institute. 
(Read at the E.R.I. March 28, 1961.—Received March 31, 1961.) 


1. Introduction 


Several years ago the same authors (MATUMOTO and SATO: 1954) 
theoretically studied the torsional oscillation period of the earth having 
a homogeneous mantle and core. Because of the assumption adopted 
that the core was a perfectly liquid or rigid body, it was presumed that 
the condition at the boundary surface was either free from stress or 
vanishing displacement. On the other hand, however, we have reports 
that the core is neither perfect liquid nor rigid (JEFFREYS: 1959; PRESS: 
1956 ; SCHLANGER: 1959), and it is desirable to know the effect of finite 
rigidity on the period and mode of vibration, so that we may examine 
the possibility of detecting the core rigidity by means of the free oscil- 
lation period. The assumption of a homogeneous mantle and core is not 
a satisfactory one in view of our present knowledge of the earth’s 
structure, but a simplified assumption will make the general nature 
clear and give us a clue for further study. 


2. Characteristic equation 


The displacement and stress components for the torsional vibration 
were given before. (Subscript o refers to the mantle and 7 the core.) 


u,—=90 
—im EV/(B I naa ir): 1 P™(cos 6)-exp (im¢)-exp (tpt 
ge n(n +1) Vr & Yasia(kor) sin 0 oe a eee 
af 1 (B drsyplkor), & —_P™(cos 7)-exp (wm¢)-ex 
p (vpt 
Dea n(n+1) Vr is we Vasa ko r) dé es ‘ “a vie 


for the mantle 
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b<r<a (2.1) 


a and b being the outer and inner radius of the mantle respectively, 
and 


0 
—im i if re , F 
= -——EJ, ..,(k,r)-———P™(cos 0) exp (im¢)- exp (7pt) 
U; nOeet Uy Fava iY) ata ( ) p ( ¢) 
w= — 4. BT, yal eer) P(cos 8)-exp (img): exp (pt) 
n(n+1) Vr dd 
for the core 
0<r<b (72) 


u, v and w are the displacement components for the radial, colatitudinal 
and azimuthal directions, and 


ky =p" Pol fo=(p/Veo) and k?=p*e,/u=(v/ Vi)? (2.3) 
where p, » and V, are the density, rigidity and S wave velocity re- 
spectively. 


The characteristic equation is given by the continuity of the dis- 
placement and stress at the boundary surface as follows 


3 0 . £-3/2 3 0 £). 6-3/2 
E FY (Snr €) 5 ) g ge etal) g ) 0 
1 — ) as 0 0 
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where 
=k, 7=k,0. and C=k,b=(Vool Vai)7 (2,4); 


3. Solution of the equation 


The above equation gives a relation between 6, b/a, k,/k, and p/p. 
For finding the effect of a thin crust over a homogeneous mantle we 
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are to assume b/a nearly equal to 1, while for the oscillation of the 
whole earth having a mantle and core we assume this value b/a=0.5447 
(=8470/6370). Then the ratio of the rigidity y,/~, and that of the 
density p,/e, are given as a function of € (=ka=2z7a/(Vol'), T=oscillation 
period) and k,/k, (= Voo/V.,). 

Fig. 1 is the result of numerical computation for f/f as a function 
€ (=k,a) and V,,/V,,, » being assumed to be 1 and 2, both axes are 
written in logarithmic scale. There are domains filled by numerals, which 
indicate the existence of solutions. On the other hand there is no 
solution for any combination of k,a and V../V., corresponding to the 
blank area, because the ratio y;,/~ (and also f:/P.) turns out to be nega- 
tive in this area. 

Horizontal lines in the figures give solutions either for a liquid or 
rigid core. L,(&,7)=0 gives the former, while R,(£,7)=0 the latter. 
These equations do not contain any physical quantity connected with 
the core, therefore the solution has nothing to do with the physical 
nature of it. A part of the result was given in the previous paper. 

Using Figs. la, 1b and two other similar figures we can draw the 
lines of constant (;/4,) and those of constant (p,/p,). They are given 
in the Figs. 2a—-2d in linear scale. (A few numerical errors were found 
in the previous paper, which are corrected here.) 


4. Discussion 


Since we can easily get the solution to the two extremity cases, 
namely the perfectly liquid and rigid cores, we are apt to think that 
the period of oscillation for finite rigidity lies between these two values. 
But the result of computation turns out not to be so. As the figures 
in the previous section show, if we keep (¢,/¢) finite 


(Oscillation period)—~, as mod. 
(Oscillation period)=finite, when 4,=0. 


If there is no core substance at all or the core rigidity is a perfect 
zero, the energy is not transmitted from the mantle into the core and 
the simple theory in the previous paper holds. However, if there is 
an interaction between them, the core has two effects to the oscil- 
lation ; namely it gives a restitutive force and also it acts as an inertia. 
The former shortens the period, while the latter elongates it. If the 
density p, remains constant and the rigidity p, tends to zero, the effect 
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of the core as an inertia term remains unchanged while the restitutive 
force becomes weak, and as a result the period becomes extremely long. 
This can be interpreted in another way. If we look at the curves 


Pi Pom 
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Fig. 3a. Curves of constant 9;/0). Bold lines refer to ¢/-axis. 
Broken lines refer to ¢-axis. 
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in Figs. 2a-2d belonging to the fundamental mode, we find (being better 
in figures of larger n) that they are divided into two parts; one part 
is nearly horizontal and is independent of V,,/V,., which fact shows that 
the frequency is approximately determined by the property of the mantle 
only. While in the other part, € is almost proportional to V,,/V,, im- 
plying that the outer part of the earth has little effect on the oscil- 
lation. Figs. 3a and 3b, which are both written using the ordinate 
&'=2ra/(V..T) (€=2za/(V,T) is also given for reference), will make clear 
the above circumstance. Whatever value », may take, the ratio of the 
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Fig. 3b. Curves of constant 0:/00. Bold lines refer to ¢/-axis. 
Broken lines refer to é-axis. 
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densities does not change widely, so the curves of constant (¢;/,) are 
given in these figures. When the shear velocity of the core is low 
curves are flat, implying the effect of the mantle being small. In this 
case the displacement in the mantle is, as is seen in the figures of the 
next section, very small, therefore the oscillation of this kind seems 
hard to be excited by a source in the upper part of the mantle. Then 
what will be the significance of the period theoretically obtained ignoring 
the rigidity of the core? The observed period of free oscillation of the 
earth agrees with this simple theory (BENIOFF, PRESS and SMITH: 1961 ; 
Ness, HARRISON and SLICHTER: 1961; ALSOP, SUTTON and EwInG: 1961). 

A similar problem was discussed by K. TAZIME (TAZIME: 1959) in 
his study on the propagation of surface waves in a medium having a 
very low velocity layer. According to what he suggests the solution 
which is obtained neglecting the existence of low velocity medium corre- 
sponds to the higher mode oscillation with such low velocity medium. 
When the core rigidity is small, higher mode has generally very large 
amplitude in the part of core and it seems that a source of disturbance 
near the surface is not a favorable condition for exciting such a mode. 
However, when the impedance contrast of two media is large the rate 
of energy transmission across the boundary is small. Besides the resti- 
tutive force at this surface is also small when the core is much less 
rigid than the mantle. Under such a condition the existence of the 
core will scarcely affect the oscillation period. Thus a mode of oscillation 
will exist, which is a higher mode if fully developed, but, in the part 
of the mantle, has at first a similar feature to the case of perfect 
liquid core. 

In Fig. 4 the left figure represents the fundamental mode where 
there is a perfect liquid core, while the right figure gives an example 
of one of the higher modes when the core has small rigidity. ,/p, and 
/4/M> Will not be very far from those of the real earth. A value of 
period (=1/é) is adopted which is nearly equal to the case of perfect 
liquid. The distribution of displacement in the mantle thus obtained 


does not differ appreciably, although the feature in the core is completely 
different. 


5. Conclusion 


The above discussion leads us to the conclusion that 


1) Existence of the core with large rigidity makes the frequency 
of torsional oscillation higher. 
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Fig. 4. Amplitude distribution 
Left: Core of liquid 
Right: Core of small rigidity 


2) Existence of the core with small rigidity decreases the frequency 
of fundamental mode and the period becomes very long. 

3) However, when the core rigidity is small, there are higher modes 
whose periods are nearly equal to those in the case of the perfect liquid 
core, and the feature of motion in the mantle is also similar. 

4) Although the displacement in the core is very large and is 
completly different from that in the case of perfect liquid, the observed 
period will be the one for these modes with undeveloped displacement 


distribution in the core. 


6. Distribution of amplitude 


To make the feature of oscillation clearer, the distribution of am- 


Yasuo SATO and Tosimatu MATUMOTO 


14 


(peeoz‘ooe) © 
(90v'0z'00z) @ 
(42v°02 ‘001) © 


(9£1z1 ‘ooe) © 
(9z0z1‘00z) © 
(60611 001) © 


LO=5 


(zs101‘00e) © | 
(€21'01‘00%) © 
(76101 001) © 


(86g ‘oor) © 
(v96's‘00z) @ 
(8z6's ‘001) @ 


vies 


(929 ‘oor) © 
(zev'9 ‘00z) @ 
(6e¥9 ‘o0'l) © 


(vize‘ooe) © 
(9zZ2 ‘ooz) @ 
(e2z¢ ‘o01) © 


(9297008) © 
(e99%‘002) @ 
(869¢ ‘oo1) © 


(vzg2‘o0g) © 
(629% ‘00z) @ 
(2692 ‘001) © 


(agee‘008) © 
(ezv'2'00z) @ 
(9ev2‘001) © 


(62s ‘o0e) © 
(8621 ‘00z) @ 
(og61 ‘001) © 


Oy= 5 


(vzgz‘o0") © 
(4z9Z‘00z) @ 
(9€2Z‘001) @ 


(ogi00‘o0e) © 
(s960 ‘o0z) © 


(o¢¢'1 001) © 


os=5 - 


n=l, 


core for 


Amplitude distribution in the mantle and the 
2nd higher mode 


Above: 
First and second numbers in the parenthes 


Big mods 


Below: 


Ist higher mode 


€S; Pi/Oo and Veo/ Vey. 


Vibration of an Elastic Globe with a Homogeneous Mantle 15 


plitude within the earth was computed for the case n=1 and 2 adding 
to the one given before. Since the density ratio does not vary widely, 
computation was carried out along the curves of P;/)=Cconst. 

As can be expected from Figs. 3a and 3b, corresponding to large 
values of V,,/V,, the energy concentration in the mantle is revealed, 


while small V,,/V,, gives a mode of oscillation which has large displace- 
ment in the core. 
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Fig. 5b. Amplitude distribution in the mantle and the core for »=2 fundamental mode. 
First and second numbers in the parentheses; 9;/09 and Vso/Vsi . 
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ABSTRACT 


1. The siress-sirain relations of some igneous rocks for longitudinal 
and lateral direction were observed. The curves of the longitudinal strain 
versus siress were approximaiely linear up io the moment of rupiure, but the 
sirain in the lateral direction showed ihe abrupi increase in the fraciure range. 

2. li was found that such abrupt increase of the lateral sirain was 
closely related with the flow. The longitudinal creep during the moderaie 
period was almost recovered by ihe removal of the load, bui in the lateral 
direction, large amount of residual sirain was observed afier the removal of 


the load. 


3. The empirical formulae of creep for granite were given as follows, 


S=Ay+ Are *"*+ Are +4 Ase “+ Blog t+Ct, 
for the longitudinal direciion, and 
S=A+ Blog t+Ct, 
for ihe lateral direction, respeciively. In ihe former equation, the recipro- 
cals of ai, a2 and as denoie the reiardaiion iimes. The laier equation con- 


curs with ihe formula given by Griggs and others. 


1. If we consider the problems of the occurrence of the earthquake, ii 
seems io be very imporiant io study che rupture mechanism of rocks con- 
stituting the earth’s crust in ihe deep place. However, ihe rocks which 
exist in the deeper pari of the crust must be different from the ones found 
at the surface of the earth. Moreover, ihey are under high pressure and 
high temperaiure. 

In spite of ihe above circumsiances, it may give us some clue of the 
inierpreiation of ihe rupiure phenomena in ihe crusi io siudy the mechanism 


of fracture of rocks in ordinary states. The experimenis reported here 


were carried out under an atmospheric pressure and room iemperaiure. In 
these experimenis ihe simple compressional force was applied io ihe speci- 
mens. 


type. The sirain gauge of this type has high sensitivity and considerable 
accuracy. Besides, it has a very convenient characierisiics such as the any 
desired strain component can be measured easily. The difficuliies caused 
by the unsiability of the indication; which we always encounter in ihe ireai- 
meni of the eleciric insirumenis during ihe long period observation, were 
avoided sufficienily by ihe use of the unbonded iype of sirain gauge. 

2. Ti has been reporied in many papers ihai ihe siress-sirain relations 
in the loading direction for silicate rocks were approximaiely linear up io 
the moment of rupiure. Ii may be said ihat the rocks are the quiie elasiic 
and briiile material. Ti is well known, however, thai ihe rocks are ihe very 
porous maierial, and because of its porous siruciure, they have ihe con- 
siderably lower value of Poisson’s raiio ai ihe incipience of loading than 
that of the other maierials?. Then, it may be expecied ai the rupiure 
siage, thai ihe porous siruciure will considerably affect io the rupture mech- 


anism for the rocks. 
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Bien 1 Stress-strain curves and variation of Poisson’s ratio with stress, (a) 
Aokiko-Nagano quartz monzonite, (b) Yakuno-Kyoto olivine basalt. 


According io the above consideration, the measurements of stress-strain 
relations of the several soris of igneous rocks both in the longitudinal and 
in the lateral directions were carried out up io the rupiure point, and ihe 
ratio of ihe laieral io the longitudinal strain was calculated. The rock speci- 
mens used for these experimenis are 30 mm. in diameter and 60 mm. long. 

The longitudinal sirain of each specimen varied almost linearly with 
siress up to the rupiure as above mentioned. As shown in Fig. 1, however, 
ic was found that the rate of increase of the lateral sirain increased gra- 
dually with siress, and in the fracture range the increase of ihe laieral 
sirain becomes much larger than the longitudinal one. At the first siage, 
ihe Poisson’s raiio gives the less value compared with the other maierials, 
and with the increasing load it reaches ai the ordinaly value, about 0.3. 
When ihe rock specimen begins io show the aspect of failure, the raiio 
increases beyond 0.5, and fimally it shows the exiraordinaly high value ai 
the rupiure siage. This suggesis ihe effect of volume expansion in ihe 
fracture or ithe plasiic stage, which was found by P. W. Bridgman in his 
volumeiric tesis”. It is supposed that this effect involves in iiself the 
fundamenial meaning for the rupiure mechanism of rocks. 

3. When the siress-strain relation was observed in the fracture range, 
it was clearly indicated that all the abrupt increase of lateral sirain was 
not produced insianianeously ait the momeni of loading, bui the flow strain 
afier loading had the close conneciion with such abrupt increase of sirain. 
Then it may be interesting io observe the amouni of flow along the laieral 
direction compared with ihe one along the longitudinal direction during the 
moderaie period at ihe fraciure stage. 

To avoid the reasonable suspicion on the barrel shape deformation of 
the specimen, the lateral strain was measured ai the ending poriion of them. 
The other conditions were kepi the same as ihe previous experiments. As an 
example, the resulis of sirain-time relations for a Kiiashirakawa biotite 
graniie under various stresses are shown in Fig. 2, where (a) denotes ihe 
longitudinal direction and (b) the lateral direction. The resulis of ihe 
others were almost the same as the above. The normal sirengih of this 
granite is 1400kg/cm?, then it will be seen that the lateral flow is hardly 
perceptible below the stress of one third of the normal strength. Beyond 
the stress of a half of the normal strength, however, the lateral flow be- 


comes larger than that of the longitudinal one, and it surpasses the later 
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Fig. 2. Strain-time curves for Kitashirakawa biotite granite under various loads, 
(a), longitudinal direction, (b) lateral direction. 


by far in the fracture range, in which ihe microscopic failure will rapidly 
occur. 

Thus the deformaiion process of rocks can be divided inio three siages, 
that is, the first stage in which the pores are closed, the second siage in 
which ihe maierial behaves almosi perfecily elastic, and the final and frac- 


ture siage in which ihe lateral sirain shows a disiinciive increase. 
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2 1000 
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Fig. 3. Stress-strain relations, which involve the flow and recovery, for Kita- 


shira biotite granite, (a) longitudinal direction, (b) lateral direction, 


In these observations, the elastic recovery was measured afier each iesi. 
Fig. 3 shows the siress-sirain-iime relation given from ihe successive test 
for flow and recovery. In ihe longitudinal direction, the sirain caused by 
ihe creep almost returned, bui the laieral creep sirain scarcely returned ai 
ihe removal of the load. Then the large amouni of permanent sei, the 
amouni of which will be closely related with ihe duration of loading, was ob- 
served in the lateral direction. Ji may be said that the rocks have the unu- 
sual properiies, elasiic in the direciion of compression and plasiic in ihe 
direction transverse io ii. 

4. Many measuremenis have been reporied on ihe creep of rocks, and 
ihe experimenial equation of ihe resulis has been given in the form of 


S=A+ Blog t+Ct, 


where S denotes the flow sirain, ¢ the time and A, B, C the consianis. 

In this section, the experimenial formula of the creep for granite is 
esiimaied. The specimens used in this experimenis are Shodoshima bioitie 
granite, 30 mm. in diameter and 60mm. long. The measurements were 
done in iwo stages, the one was in ihe short period 1 io 10 second, and 
the other was in ihe long period 10 io 10° seconds. The creep of long 
period was measured by the unbonded-iype sirain gauge. _ This sensiivity is 
a quariver of 10-° in the axial direction and 10-° in the transverse direciion. 
The error of the observed value was within 3x 10-5, during the all period 


of ihe observation. The change of ihe room iemperaiure, by which the 
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Fig. 4. Creep curves of longitudinal direction for biotite granite, 

test for Kitashirakawa granite under load of 945 kg/cem2, (b 
for Shodoshima granite under load of 1170 kg/em?. 
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observed values are maierially affected, was kepi within +0.2°C in the axial 
direction and +0.5°C in the transverse direction. The resulis are described 
separaiely in two direciions. 

In Fig. 4(a) and (b) show the creep curves for the longitudinal direc- 
tion in two periods ploiied on semilogarithmic papers, respeciively. The 
empirical formulae derived from these iwo curves are 

S=ALA'le “Ci, 

S=A)t+Aie “!+4 Are 4 Blog t+ Ct. 
Let us assume that ihe term Cf? of ihe former equaiion can be replaced by 
the second and laier terms of the laier equaiion, as it was deiermined in 
the short period | io 10? seconds. Then the full equation can be described 
in the form 

S=AotAie “44 Ase” 44 Age + Blog t+Ct 

puiiing together, where the reciprocals of a1, G2, as express the reiardaition 
times of Kelvin’s model of the order of 10, 10°, 10‘ seconds respectively. 


The resulis of laieral direction in shori and long periods are shown in 
Fig. 5. These formulae are 


ge ee THe? io! 
(a) Dime (sec) (b) 
Fig. 5. Creep curves of lateral direction for Shodoshima biotite granite under 
load of 1170 kg/cm?, (a) short period test, (b) long period test. 


S=A+Blogt, 
S=A+Blog t+ Ct, 


respeciively. Putting iogeiher, ihe empirical formula of ihe creep for the 


8 


direction perpendicular io the applied force is given as follows 


S=A+Blogt+Ct. 


This formula entirely coincides wiih ithe formula given by Griggs and 
ozhers”. The formula for ihe longitudinal direciion, does noi concur with 
theirs, but has ihe different expression with the lateral one. While ihe 
curves shown in Fig. 5 are approximately linear, the ones in Fig. 4 some- 


what indicate the unduraiion. 


, Tupture 


\ appreciable crack 


Strath 


0 10 20 30 40 50 60 70 
Time (hour) 
Fig. 6. Creep and rupture for Shodoshima biotite granite, observed in lateral 
direction. Solid line denotes the observed value, and the dotted line denotes 
the normal creep curve. 


An example of rupiure caused by creep is shown in Fig. 6. This curve has 
the same three siages as the meials and other materials have. However, it 
was observed that the deformation did noi continue smooihly, bui the ap- 
preciable cracks occurred in succession on ihe way of creep, and finally 


reached io rupiure with continual crack. 
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Abstract 


It seems that the strength of the earth’s crust and upper mantle has 
some relation with the magnitude of the earthquakes. In this paper, we 
try to estimate the strength distribution of the earth’s crust and upper 
mantle from the experimental results carried out under high pressure and 
tem perature. 

We compare the above result with the distribution of the great ear- 
thquakes with depth, which occurred in and near Japan in 1926-1956. It 


seems that the both are fairly resemble each other in tendency. 
Introduction 


The geographycal distribution of the earthquakes has been studied 
in detail and well known as the seismic zone. On the other hand, notic- 
ing the distribution with depth, it will be seen that the distribution of 
the hypocentres of the great earthquakes with depth has the considerably 
conspicious feature, even of the shallow focus earthquakes within the 
depth of 100 k.m.. 

The abundant problems pertaining to the mechanism of the earthqua- 
kes will not be able to be proved till we get the sufficient knowledges for 


the time and space distribution of the forces in the crust and upper 


mantle, and for the constitution of there. Nevertheless, it will be the 


reasonable inference that the magnitude of the earthquakes will depend 
on the mechanical strength of the hypocentre, and it will not be the 


futile attempt to estimate the strength of the crust and upper mantle. 


Because the large amount of energy should be stored in the region where 
the strength is large, and the less do in the region of the small strength. 

In this paper, we try to estimate the distribution of the strength of 
the crust and upper mantle with depth by the experimental results on 
the strength of rocks under high pressure and temperature. Then we 
consider the distribution of the great earthquakes with depth, which oc- 


curred in and near Japan, with references to the above estimated results. 


Pressure, Temperature and Constituents 
of the Crust and Upper Mantle 


The strength of rocks is strongly affected by confining pressure and 
temperature, and these effects are different for rocks of different types. 
Then we must know the 


pressure and tempera- 


ture distributions and 1500 30 
the constituent rocks of 
the crust and upper @ g 
5 1000 20 7 
mantle to evaluate the % a 
strength there. O = 
: e 
Of the pressure dis- S00 10 


tribution, we have had 


the sufficient knowledge 
1) 20 40 60 80 100 km. 
Depth 


Fig. 1. Pressure and temperature distribution in the 


for our purposes’. This 
is shown in Fig. 1. 
enor cmec pee Oat earth’s crust and the upper mantle. 
distribution, on the other 
hand, the various estimations have been presented by many authors” and 
we have had no definitive distribution. These distribution curves, how- 
ever, are approximately coincident in tendeny and magnitude one another. 
So we will take the average of these curves as the temperature distribu- 
tion of the crust and upper mantle. This curve is shown in Fig. 1, which 
will not give so large error for the temperature within the superfical 
layers of the earth. 

Next we consider the constituent rocks of the crust and upper mantle. 


The boundary between the crust and the upper mantle shows very re- 
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markable discontinuity and has been well "known as Mohorovicie dis- 
continuity. From the seismological and geological evidences, it has been 
clarified that the crust is composed of relatively acidic rocks such as 
granite and that the upper mantle is done of basic rocks such as dunite 
or so. In the following discussions, we assume the simplest model of 
the crust and upper mantle, that is, the crust is composed of granite and 
the main constituent rocks of the upper mantle are dunite or similar 


rocks. 


Variation of the Strength of Rocks with Confining Pressure 


In the previous paper”, we studied the variation of strength with 


various confining  pres- 
% sures for Kitashirakawa 
granite. This and subse- 
20 quently obtained results 


are shown in Fig. 2. The 


empirical formula of the 


Strength 


strength versus pressure 


(=) 


for this rock is given 
p* = Po* (Re p+1) 1, 
where Po¥, Pp* are 


strength at an atmospheric 


R'=15 (1722p +1)4 kb. 
the 


0 5 10 


15 kb. 


pressure and 


Confining Pressure 


Fig. 2. Variation of strength-of granite with 
confining pressure at room temperature, 


confining 
pressure’ of 1p. k. kpere- 


spectively, and & is the 


characteristic constant for 
the rock. 


Here we assume that the above formula holds in all types of silicate 
rocks, though we can. not assert strongly this assumption for the accumu- 
lation of the experimental results are not sufficient to determine the 
plausible formula. 

As k and P)* for granite, we use the experimentally derived values 
for Kitashirakawa granite, that is, R=17.22 (k.b.)-! and Py*=1.5 kb. 
The average value of P*, listed in ‘‘Handbook of Physical Constants®”’, 


is 1.48 k.b., so that the above value seems to be appropriate. 


Next, we must find the values of k and P* for dunitic rocks, There 
has been no experimentally derived value of k for this sort of rocks, so 
we will determine this value from the values of Pp* and P,* (strength 
under certain confining pressure). In “Handbook of Physical Constants’’, 
average value of P)* for gabbro, diabase, etc. is listed. This is 1.8k.b.. 
The strength of gabbro, basalt and diabase, which we obtained, was about 
2.0 k.b.. Then we put Po*=2.0 k.b. for dunific rocks as the approximate 
value. D. T. Griggs et al.’ has observed the stress: strain relations for 
dunite under 5.0 k.b. confining pressure. In their report, the strength 
of dunite at room temperature is shown as about 21 k.b.. From the 
values of P,* and P;*, we can obtain the value of k=22.0 (k.b.)- for 
dunitic rocks. 

In Table 1, the strength of granitic and dunitic rocks under pressures 
corresponding to the various depths from the earth’s surface, which is 


calculated by the above mentioned way, is listed. 


Table 1. Strength of granite and dunitic rocks under pressures corresponding 
to the depths from the earth’s surface. 


Strength 
Depth Pressure é 
Granite Dunitic rocks 
km. oe Keb Ps. oa k.b. 
0 0 is 2.0 
5 13 (or 10.9 
10 Ath 10.3 15,5 
20 5.2 14.7 22.1 
30 8.5 18.2 27.4 
40 is 74 Wis ay 
50 15.0 24,2 36.4 
60 18.5 26.8 40.4 
80 Ae ice 47,4 
100 32.0 See, 53. 1 


Reduction of the Strength caused by the 
Elevated Temperature 


It has been currently accepted that the strength of rocks is consider- 
ably reduced by the elevated temperature, but it has not been made 


clear quantitatively how much amount of reduction is brought about. 


Recently, D.~ T. 


Griggs et al., as above 


kb. 


mentioned, has  experi- 50 


mentally investigated the 


stress-strain relations for 


& 
various sorts of rocks at %!5 
S 
the temperatures up to @ 
ce 
800°C and pressures of H a 
(0) Veloye 
We will quote their 
results and picture again 5 
the temperature-strength 
relations for granite and 0 
dunite. This is shown in 0 
Fig. 3. As seen in this } 
Pigase 


figure, it seems that the 
linear relation exists bet- 
ween the strength and tempera- 
ture for granite, and exponen- 
tially decreasing relation for 
dunite. Here we assume that 
the following formula can be 
adopted for granite as the em- 
pirical one and that this formula 
is held all over the pressure 


range, that is, 


pPy* = >Pp* (1 = T/1100) Py 


Fig. 4. Variation of strength of 
granite and olivine basalt with 
temperature, at an atmospheric 
pressure. The upper and lower solid 
curves show the crushing strength 
of olivine basalt and granite re- 
spectively, and the dotted curve 
shows the estimated creep limit of 
granite. 


Granite 


Dunite 


400 600 800 1000 
Temperature 


1200°C 


Variation of strength of granite and dunite 
with temperature, at 5.0 k.b confining pressure 
(after D. T. Griggs et al.). 


Olivine Basalt 
Crushing Strength 


a 


Granite 
# Crushing Strength 


NX 
5 \ 
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Temp erature 


where 7P,*, >Pp* are the strength at temperatures of T°C and 0°C under 
confining pressure p k.b., respectively, 

The temperature-strength relation for dunite is not given in such a 
simple form as for granite, then we evaluate the strength reduction 
graphycally from the curve in Fig. 3. That is, the strength at T°C, p 
k.b. is given as the product of the strength at O°C, p k.b., which was 
estimated in the former section, by the ratio of strength at T°C and 0°C 
under 5 k.b.. This is given formally as follows, 

rPy* = Py* (7P3*/oPs*) : 


In addition to the above argument, we refer to the temperature de- 


Strength 


0 20 40 60 80 100 k.m. 
Depth 


Fig. 5. Strength distribution in the earth’s crust and the upper mantle. 


Table 2. Strength of the earth’s crust and the upper mantle taking account of 


the temperature and pressure. 


| Strength 
Depth Pressure Temperature | 7 es 
Granite | Dunitic rocks 
eee ey Kb {coe k.b | 
| 0 0 1.5 2.0 
5 1,3 200 5.95 8.0 

10 Pi 350 7.05 9.15 | 

20 | 5.5 550 | 7,35 | 10.15 
| 30 8.5 700 6.6 10.6 

40 11.5 850 4.8 10,15 

50 15.0 980 | 2.6 9.9 

60 18.5 1100 | 0.0 9.55 

80 25.5 1300 | 8.2 
| 100 | 32.0 1500 6.3 
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pendence of strength of granite and olivine basalt in our experiments at 
an atmospheric pressure®. This is shown in Fig. 4. This results will 
suggest that the above treatment is not so unreasonable. 

Thus we can estimate the strength distribution within the earth’s 
crust and upper mantle taking account of the temperature and pressure 


of there. This is shown in Table 2 and Fig. 5. 
Creep Limit of the Earth’s Crust 


If the earthquakes are caused by forces generated suddenly®’, the 
crushing (or yield) strength obtained in the foregoing sections should 
give the actual one in the crust and upper mantle. On the other hand, 
if the earthquakes are caused by forces accumulated gradually, the above 
obtained strength does not give the real strength. The stress is released 
constantly by the flow or the creep in such a rather static process, and 
rupture must be occurred by further lower stress than the crushing 
strength. Then we must estimate the creep limit in the crust and upper 
mantle. 

For this problem, we scarcely have any experimental data to be avail- 
able. Nevertheless, we will try to estimate the rather qualititative dis- 
tribution of the creep limit within the earth’s crust, using the very 


rough-and-ready presump- 


tions. Creep limit can be 
obtained only by the very 
long period tests, but such 
experiments are considera- 
bly difficult to be carried 


out sufficiently, especially 


at high pressure and tem- 


Rate of Creep(Arbitary Scale) 


: o perature. Thus we assume 
reep Limit 


that the creep limit cor- 


0 500 1000 ae . responds to the point at 

Stra Kg clans 3 
alia which the curve of the 
Fig. 6. Relation between the rate of creep and 


stress, and the estimated creep limit for Kita- 
shirakawa granite, at an atmospheric pressure comes In Contact asym pto- 


and room temperature, tically with the stress axis. 


rate of creep versus stress 


In Fig. 6, the curve of the rate 
of creep versus stress for Kita- 
shirakawa granite at an atmos- 
pheric pressure and room tem- 
perature is shown. Fig. 7 shows 
the variation of creep limit with 
Kita- 
shirakawa granite obtained by 
The dotted 


line in Fig. 4 shows the relation 


confining pressures for 


the above method. 
between creep limit and tem- 


perature at an atmospheric pres- 


sure for this granite. 


Creep Limit 


0 1 2 3 4 5 Gaukibe 
Confining Pressure 


Fig. 7. Variation of the estimated creep limit 
with confining pressure for Kitashrakawa 


granite, at room temperature. 


It will be seen from the figures that the empirical formulas of the 


variation of creep limit with confining pressure and temperature have 


the same form as the one of crushing (or yield) strength. 


ivateniss 


Pee hae(k-p+ la, 


and 


il — wet (1 = T/800) 9 


where the index c instead of * appeared in crushing (or yield) strength 


denotes the creep limit. 

Using the results shown in 
Fig. 7 and Fig. 4, we can rough- 
ly estimate the creep limit of 
the earth’s crust. The result 
is shown in Fig. 8. Of course, 
the 


result must be recomputed after 


quantitatively plausible 


the accurate tests of creep under 


high pressure and temperature. 


kb 

2 
i 
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O 
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Fig. 8. Distribution of the creep limit with 
depth in the crust. 


Distribution of the Great Earthquakes with Depth 


Let us consider the distribution of hypocentres of the great earth- 


quakes with depth. 


We pick up the shallow focus earthquakes, the 


magnitudes of which are larger than 6.3, which occurred in and near 


10 


Japan from 1926 to 1956. The data are listed in “The Seismological Bul- 
letin of the Japan Meteorological Agency, Supplementary Volume Nov 
uae 

We consider the seimic zone in Japan dividing into two parts by 
Fuji Volcanic, Belt. In North-eastern Japan, the earthquakes have occur- 
red mainly under the Pacific Ocean along the coast and these hypocentres 
are rather deep in general. On the other hand, in South-western Japan, 
the earthquakes have occurred mainly under the land and the depth of 
hypocentres is shallow. This means, that is to say, that the earthquakes 


have occurred in the upper mantle in North-eastern Japan and in the 


Fig. 9. Geographycal distribution of the great earthquakes which occurred in 
and near Japan in 1926-1956. 
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Fig. 10(b). 

Fig. 10. Distribution of the great earthquakes with depth, which occurred in and 
near Japth, in 1926-1956. (a) In North-eastern Japan, and (b) in South-western 
Japan. The full circles and empty circles express th magnitudes determined by the 
Japan Meteor. Agency and by B. Gutenberg respectively, and the X marks express 
the earthquakes, the depths of which are not determined definitely. 


Fig. 9 shows the geographycal distribu- 


crust in South-western Japan. 
Fig. 10(a), (b) show the distribu- 


tion of epicentres in and near Japan. 
tions of hypocentres in North-eastern and South-western Japan respectively. 
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Comparing Fig. 10 with Fig. 5, it seems that the tendencies of both 
distributions with depth are fairly alike to each other. Fig. 10(a) cor- 
responds to the upper curve in Fig. 5, strength of dunitic rocks, and 
Fig. 10(b) to the lower curve, strength of granite. 

If the energy of earthquake is chosen as the vertical coordinate, the 
curves in Fig. 10 show the rather sharp maximum. This will mean that 
the concept of the “earthquake volume®”’ is very available, because the 
earthquake which occurs in the region with great strength should have 
the much larger volume than the one which occurs in the region with 


less strength. 
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1. Introduction 


For many years microseisms have been investigated by numerous 
scholars. But there is still lack of general agreement about the origin, 
the mechanism of generation and the other properties of microseisms. It 
is the object of this paper to solve these problems. 

It has been well known that microseismic storms appear on seismo- 
grams during the passage of typhoons or cyclones. The view has often 
been expressed that microseisms are associated with atmospheric depres- 
sions, but recently most of investigators support the opinion that sea 
waves generate microseismic waves in the solid crust. 

Bernard (1937) et al. found that the period of microseismic waves is 
about half that of the generating sea waves. Longuet-Higgins (1950) 
presented a significant theory for the effect at the ocean bottom of stand- 
ing sea waves produced by the interference between two similar wave 
trains travelling in opposit directions. Such a condition as Longuet- 
Higgins pointed out is expected to be attained at the neighbourhood of 
the atmospheric low pressure center or the coast, and hence the opinion 
on the location of the origin of microseisms is divided into two.  Al- 
though at present his theory seems most attractive, no evidence of ex- 
istence of standing sea waves is observed. 

From the case studies of microseismic storms, it is known that the 
occurrence of the peak in microseismic activity observed at an obsery- 
ing station does not coincide with the moment when a typhoon or cyclon 
approaches to the station as closely as possible, but rather is delayed, 
and also the occurrence of the maximum amplitude at each station seems 


to propagate with the velocity similar to that of progress of the typhoon 
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or cyclone. To explain the phenomenon some of seismologists considered 
that the generation of microseisms is due to the swells propagated from 
a storm area to the coasts, and they investigated the relation between 
microseisms and swells. Sakata (1940) presented an excellent paper on 
this subject. Wadati et al. also were led to the same consideration from 
the existence of close relation between microseisms and swells, as observ- 
ed at some stations in Japan. Santo (1959) also investigated in detail 
the problem by the observational data obtained at many stations in Japan 
during I.G.Y. period and confirmed that consideration. 

One of the most effective means to investigate the origin of micro- 
seisms is a determination of arrival directions of microseismic waves. In 
spite of a great deal of efforts of seismologists, they did not reach the 
satisfactory conclusion to that investigation because of their unsuitable 
methods, Therefore the writer did not follow their methods, namely, 
the tripartite and the Lee’s method, and he studied the arrival directions 
from the analysis of orbital motions of the earth’s particles by means of 
vector seismographs. It became evident from these observations that the 
arrival directions of microseismic waves are associated with not the posi- 
tions of atmospheric low pressure centers, but the coasts near the ob- 
serving station. And a comprehension of the other properties of micro- 
seisms was also considerably gained, and these results appear to support 
the above-mentioned consideration as already pointed by the writer (1959) 


and (1960). 


2. General features of microseismic waves observed 
at the Abuyama Seismological Observatory 


The writer will at first refer to the general features of microseismic 
waves at the Observatory based upon the information deduced from the 
data on the waves by the precise observation during I.G.Y. period. The 
Abuyama Seismological Observatory is located at the position with co- 
ordinates : 34°52’N and 135°34’E, and founded on palaeozoic system. 

Good’care was taken in making selection of the constants of seismo- 
graphs. With due consideration for ground movements with periods of 
2 to 8 sec. and with amplitudes up to 7 microns, the constants of seismo- 


graphs were arranged as follows. 


Constants of seismographs. 


Component | ie as hy hy Vmax. 
UD 4.00 10.00 0.%5 1.00 4,100 
NS 4,00 10,00 OFS 1.00 3,500 
EW 4,00 | 10,00 Onto 1.00 3,400 

Lan : Period of pendulum 

T2 : Period of galvanometer 

hy : Damping constant of pendulum 

he : Damping constant of galvanometer 


Vmax.: Maximum magnification 


Wave type 


From most observations hitherto made in the world it was reported 
that the three components of microseismic waves possess nearly same 
amplitudes, and hence it has been widely accepted that microseismic waves 
are the Rayleigh waves. Some of seismologists, however, have considered 
them to be of the Love type, the Rayleigh type combined with the Love 
type or the Rayleigh type combined with the standing wave. The micro- 


‘a WEE ce ee 4 ( ce wo fk e Be 
Vette Oe ee a f 
cS 5 te ee Cee SOue Sel ie 


Bie GL). 


4 oh Nei \ 7 re - \ oN 


Ra 149 


Ns 


146 


Uist. IGA). 


— 
> 
op) 
4 
4 
=> 
a 
> 
wees 


HOO GAN 
168 69 Pao 170 {TI ; 


Hig Gi) 


Fig. 1. Examples of seismograms. 


seismic waves observed at the Abuyama Seismological Observatory appear 
to be of the Rayleigh type, as shown in Fig. 1, which is the example of 
seismograms obtained by the vector seismographs. In fact the vector 
seismographs recorded frequently the waves of the typical Rayleigh type, 
but on the other hand, recorded little the Love or standing waves, 
apart from a few waves regarded as apparent Love waves or standing 
waves produced by superposition of several waves propagated from various 


directions. 


Amplitude 


The ground amplitude and the wave period by which the writer will in- 
vestigate microseisms in following pages were measured according to the 
Instruction Manual of Seismology for I.G.Y. 

Fig. 2 shows the frequency distributions of ground amplitudes in 
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Fig. 2. Frequency distributions of the amplitude of microseismic waves. The 


numerals show mean values of amplitudes in each component. 


three components. The numerals in the figures give the mean values of 
amplitudes in each component. The value in vertical component is similar 
to that in NS-component, whereas the value in EW-component is somewhat 
less than in NS-component. This tendency of amplitude is observable also 
on the other seismographs of the Observatory. If the microseismic waves 
are of the Rayleigh type and come from the coasts near the Observatory, 
this tendency should be natural and resanable, as there is no coast close to 
the Observatory in its easterly and westerly regions, but in the northerly 
region the Sea of Japan and in the southerly region the Pacific Ocean are 


comparatively close to the Observatory. 
Period 


Frequency distributions of wave periods are shown in Fig. joe Lhe 


numerals in the figures give the mean values of periods in each com- 
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Fig. 3. Frequency distributions of the i i i 

3. Freq s period of spective max 

microseismic waves. The numerals show mean P TUE 
values of periods in each component. frequency in the dis- 


tributions. The group whose peak is at the period of 4 sec. is what you 
call microseisms. Another group with the peak of 2 sec. has been ob- 
served at various stations during the passage of cold fronts. Lynch (1952) 
studied this 2-second microseisms observed at New York by using the 
tripartite observation. The results obtained by his investigation indicate 
that the Great Lake is probably the source of generation OF the 2-second 


microseisms. This sort of microseisms will be referred afterwards 


3. Variation of the amplitude and the period of microseismic 


waves during microseismic storms 


When typhoons or cyclones are passing near an observatory, micro- 
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Fig. 4. Ralation between amplitudes and periods of microseismic waves, 


seismic waves grow in amplitude and their period changes in the wide 
range, There are many reports of these variations, but they gave the 
explanation only for the case of the individual storm. The writer studi- 
ed these variation in detail during I.G.Y. period and could give the 
plausible explanation of them. Since microseismic waves resemble in 
appearance in three components, hereafter they will be stated only in 
NS-component. Fig. 4 shows the relation between ground amplitudes 
and wave periods of microseisms. The points are scattered, but they 
have the tendency that the amplitude is larger for a longer period 
as pointed out by most of seismologists. To investigate this problem 
more precisely than made hitherto, the writer gave attention for the 
case of the individual storms. The symboles excepting the solid 
circles in Fig. 4 show the waves whose forms are remarkably. re- 
gular in appearance in all seismogrms. They fall into three groups 
as indicated by the broken curves, that is, the group laid along the 
abscissa (A), near the middle region (B) and near the upper limit 
of the scattered points (C). Then the writer made a comparison between 
those three groups and the meteorological conditions corresponding to the 
times when the three groups were observed. Fig. 5, 6 and 7 are the 
weather charts corresponding to the times when the waves expressed by 


the symbols of open triangles (A), open circles and crosses were observ- 
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circles in Fig. 4 were observed. 


ed respectively, 

The first group (Fig. 5) was observed at the time when the typhoon was 
far distant from the Observatory, the second (Fig. 6) during its presence 
at a moderate distance and the third (Fig. 7) on its passing near the 
Observatory. From this fact the following inference may be drawn. Sea 
waves generated in a typhoon area have the wave height and period de- 


ending on the scale of typhoon, After leaving the @enerating area they 
p g yt} £ g g ) 


Fig. 7. Weather charts at the time when the waves expressed by symbols of crosses - 
in Fig. 4 were observed. 


propagate through the open sea to the coasts as the swell, and generate 
the microseismie waves in the neighbourhood of the coasts. The oceano- 
graphy suggests that the swell decreases in height and its period becomes 
long with propagating of 


the swell. Therefore the 250 km 
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Fig. 8 shows the re- Period of swell (sec.) 

Fig. 8. Distance (km) from which swell comes as 

bevetr amd the period of viaaiiiiag of height and period of swell at end of 
distance of decay. 

swells at the end of dis- 

tance of decay given by Sverdrup and Munk (1947). A similarity bet- 


ween Fig. 4 and 8 is very noticeable and this is one of the probable 


lation between the wave 


proofs that the mictoseismic waves may be generated by the swell pro- 


pagated from a typhoon area to the coasts. Accordingly the scale of ty- 
phoon (wind velocity in generating area, duration of wind and so on) 


and the decay distance of the swell, that is, the distance from the typhoon 
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to the coast near the station, may be estimated by using Fig. 4 from the 
amplitude and the period of the microseismic waves at the Observatory. 
There is another group of waves with the period of 2 sec. in Fig. 4. 
This group also has the tendency growing in amplitude with an incréase 
of the wave period and the slopes by which those two groups are envelop- 
ed resemble closely. But the group with the period of 2 sec. is lacking 
in points along the abscissa. Accordingly it is supposed, that the 2-second 
microseisms are probably generated by the sea waves that the local winds 
due to the cold front produce in the neighbourhood of coasts, while the 
4-second microseisms are generated by the swell propagated from the 
disturbance source far away from coasts. 

Apart from the general discussion of the amplitude and the period of 
microseismic waves as described above, the writer will turn to the discus- 
sion on the individual character of them referred to the two kinds of 


storms, e.g., the cyclone in winter season and the typhoon in summer 
season. 
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Fig. 9. Variations of distance from the Observatory to center of depression (A), 


central pressure of depression (B), periods (C) and amplitudes (D) of micro- 
seismic waves. 


Fig. 9 shows the variation of the amplitude and the period of micro- 
seismic waves generated by cyclones in winter season, the central pressure 
of depression and the distance from the Observatory to the center of 


depression. And Fig. 10 shows the variation of microseismic waves 
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Fig. 10. Variations of distance from the Observatory to center of typhoon (A), 
central pressure of typhoon (B), periods (C) and amplitudes (D) of microseismic 
waves. 


generated by typhoons in summer season, the central pressure of typhoon 
and the distance from the Observatory to the center of typhoon. It is 
marked by three noticeable differences between those two variations. The 
first is that the wave period in winter is shorter than in summer, that 
is to say, in winter the period variates centering around about four 
seconds, and in summer around about five seconds. The second is the 
marked growth of the amplitude against the gradual growth of the period 
in winter, while in summer the growth of the amplitude is similar to 
that of the period. The third is that the fall of growth of the amplitude is 
earlier than that of the period in winter and it is the opposite in summer. 
Since the seasonal wind stirs mainly the Sea of Japan in winter and the 
typhoon carries swells from a distance to the Pacific coasts.in summer, 
the disturbance source in winter is not so far as the occasion of the ty- 
phoon. Accordingly the wave period in winter is shorter than in summer 
when the microseismic waves are produced by the swell with a longer 
distance of propagation, and it rises not frequently above five seconds. , 
In winter season a low pressure produced in the south-westerly ocean off , 
Japan passes over the Sea of Japan with developing its energy, as the:, 
variations of the central pressure of depressions and the distance from the,{ 


Observatory to the center of depressions are indicated in Fig. 9, and hence, 
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the microseismic amplitude increases rapidly and its period increases 
moderately by the decrease of the distance of propagation of swells and 
the development of the storm energy. When the disturbance source is 
going away, the amplitude decreases by the increase of the distance of 
propagation of swells and the period increases furthermore by the de- 
velopment of the storm energy and the increase of the distance of propa- 


gation of swells. 
On the other hand the energy of a typhoon produced on the souther- 


ly Pacific Ocean far off Japan at first increases with the northward ap- 
proaching of the typhoon as shown in Fig. 10. Therefore the microseismic 
waves generated by the typhoon grow gradually in amplitude and period 
owing to the increase of the storm energy which is expected to be power- 
full enough to exceed the opposite effect due to the decrease of the distance 
of propagation of swells. ‘As the typhoon approaches considerably near 
Japan, the energy begins to decrease. Then the wave period of micro- 
seisms also begins to decrease owing to the decrease of the energy of 
typhoon and the diminishing of the distance of propagation of swells, 
wherease the amplitude rises’ still owing to the diminishing of the distance 
of propagation of swells. It has been hitherto considered that the wave 
period of microseisms should be associated with the depth of the ocean 
where a disturbance source exists. But this offers a striking contrast 
with the result, as above mentioned, that the dependency of the wave 
period on the energy scale of the disturbance source and the distance 
from the disturbance source to the coast near the observation point is 
obviously proved from the writer’s observation. As seen on the figures 
of Santo’s paper (1959), the variations of the amplitude and the period of 
microseismic waves depend a little on the location of the observing 
station, though the amplitude is influenced by the vibrational properties 
of the ground and the period by the instrumental constants of seismo- 
graphs. From this fact, it may be accepted that the observational data 
of most of stations located at least in Japan may yield the same results 
as the above one deduced from the observation of the Abuyama Seismo- 
logical Observatory alone. The wave period of about 4 sec. which pre- 
dominates in microseisms observed in most of stations of the world, has 
been frequently regarded as being due to the vibrational properties of 


the ground. However, if it is so, the fact that there is no maximum peak 
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of the amplitude corresponding to the period of 4 sec. can not be accepted 
without contradication (see Fig. 4). 


4. 


Tripartite observation 


x 


The tripartite observation also was carried out during I.G.Y. period, 


but the satisfactory result was not 


brought. The positions of three points : 
of observation are shown by A, B and ; 
C in Fig. 11. The instruments used t 


in this observation are the vertical 


electromagnetic 


were adjusted to the following con- 


stants: 


stant of pendulum are 4.0 sec. and 0.65, 

the period and the damping constant 

of galvanometer are 10.00 sec. and 1.00 

and the magnification is about 4,000. N 
As the side-lengths of the tripartite 101m 


net were not sufficiently long notwith- 


seismographs 


the period and the damping con- 


which Seong 


860m 


Positions of observation. 


Fig. ii. 


standing the serious attention for the 


coincidence of constants 
of the 
graphs, the considerable 


three seismo- 
error was unavoidable. 
The velocities of propa- 
gation and the arrival 
directions gained by the 
observation were scatter- 
ed in the wide range. 
The relation between ve- 
locities and periods is 


shown in Fig. 12. About 
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Fig. 12. Relation between velocities and periods 


of microseismic waves. 


half of the numbers of velocities computed are above 3 km/sec, and they 


are excluded from Fig. 12 because they are of unreliable accuracy. The 


open circles show the values of the waves which have regular sinusoidal 
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Fig. 13. Arrival directions of microseismic waves measured by tripartite method. 


forms and whose periods are identical within 0.1 seonds at three points 
of observation. The arrival directions of individual waves are shown in 
Fig. 13. The directions are considerably scattered, but it may be seen 
that they point toward mainly the Sea of Japan and the Pacific Ocean. 
And we notice the tendency that the velocity increases with the decrease 
of the period. Ikegami and Kishinouye (1951) also indicated this tendency 
in their paper, but the subject does not bear further discussions because 
of the insufficient accuracy of the observation. As Don Leet (1949) dis- 
cusses, Many microseismic waves coming from various directions are su- 
perposed and the each phase of the superposed waves differs at different 
points of observation. Therefore it stands to reason, that the values of 
velocities and arrival directions of microseismic waves are scattered in 
wide range, if they are computed from the phases selected at will in the 
tripartite method. The writer can hardly avoid the conclusion from his 
observation that the tripartite measurement is not effective to the de- 


termination of the velocity and the arrival direction of microseismic waves. 
5. Observation by vector seismographs 


From the above mentioned observations it is a natural inference that 


the microseismic waves are generated by swells propagated from the dis- 
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turbance source to the coasts near the observing station. But this in- 
ference was not drawn from the direct observation of the origin of 
microseisms. It seems most important to the writer that the arrival 
direction of microseismic waves is clearly observed, if the best results 
are to be obtained for the investigation of the origin of microseisms. 
For this purpose, the writer adopted the analysis of the orbital motions 
of waves. The orbital motions analyzed hitherto by many seismologists 
were found to be so complicated that any usefull information might not 
be deduced. To overcome this difficulty the writer made newly vector 


seismographs. The seismograms re- 


; ; Mirror of 
corded at the Abuyama Seismological galvanometer 


Observatory were also considerably Prism 
complicated, but their orbital motions 
brought out the various significant 
facts. The recording system of orbital 
motions is shown Fig. 14. 

The seismographs using for the Light 
routine observation during I.G.Y. Seen 
period were turned to immediate ac- 


count for this purpose, and hence the Fig. 14. Recording system of orbital 


constants of seismographs are the same jae 

as on the occasion of the routine observation. Three orbital motions in 
UD-EW, UD-NS and EW-NS planes were simultaneousely recorded on 
the same recording paper which runs during 1 sec. and stops during about 
4 sec. Examples of the seismograms are shown in Fig. 1. The single 
wave of the Rayleigh type is frequently observed, but the Love wave or 
the standing wave is scarcely present. Only a few of the later two wave 
types may be considered to have been accidentally produced by the super- 
position of waves coming from various directions. On comparing the 
seismograms of ordinary seismgraphs with those obtained by vector seis- 
mographs, it may be readily understood that even if the vertical amplitude 
is small and the horizontal amplitude is considerable or even if they are 
the opposite on the ordinary seismographs, it is in consequence of the 


ae i i ismic waves are regarded 
superposition of waves. Accordingly the microseismic wa g 


as the Rayleigh waves in the following discussions. 
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Frequency distributions of arrival directions 


The longer axis of the elliptical orbits of ground particles in the 


horizontal plane was used for the study on the direction of approach of 


microseismic waves by Donn (1954), Strobach (1955) and others, but their 


inferences are insufficient for lack of the records in the vertical com- 


ponent. The writer picked out such waves that their orbits are nearly 
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linear in the horizontal plane 
and elliptic in the two vertical 
planes for making the frequen- 
cy distributions of arrival di- 
rections. ©The distributions 
were directly drawn on the 
charts which were divided 
radially into thirty-six round 
the Observatory. Fig. 15 
shows the distribution on the 
occasion of the microseisms 
generated by the typhoon No. 
6 in 1959 and the travelling 
path of the typhoon is shown 


in Fig. 16, 
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Fig. 15. Frequency distributions of arrival directions of microseismic waves. 


When the center of 
the typhoon was passing 
off the south-west coast of ve 
Kyushu the microseismic 
waves come to the Ob- 
servatory mainly from 
the Pacific coasts as 
shown in (a), (b) and (c) 
of Fig. 15, and their 
directions of propagation 
do not point toward the 
center of the typhoon. As 
the typhoon did not yet 


stir up the Sea of Japan, 


the waves appearing to 


come from the coasts of Aug. 
cdnoh (7 ¢ a 
the Sea of Japan are Fe CESOR! CZ. Ser 
scarcely found. When the Fig. 16. Travelling path of the typhoon No. 6 


typhoon passed through in 1959. 


Kyushu and went on off the east coast of Kyushu, the distribution is (d) 
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of Fig. 15. It is similar in appearance to (a), (b) and (c), and the ar- 
rival directions do not point toward the center of the typhoon. When the 
center of the typhoon went on into the Pacific Ocean on Aug. 10, the distri- 
bution is shown in (e) of Fig. 15. In this case the Sea of Japan was 
stirred up and in the southern ocean off the Kii Peninsula swells went 
down considerably, and hence most of the arrival directions point toward 
the Sea of Japan. In winter season the seasonal winds blow hardly over 
the Sea of Japan and the surge is raised, whereas the Pacific Ocean 
remains comparatively calm. Therefore the frequency distributions of 


the propagating direction of microseismic waves pointed mostly toward 
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Fig. 17. Frequency distributions of arrival directions of microseismic waves. 


the coasts of the Sea of Japan as shown in Fig. 17, When the disturbance 
sources are in the Sea of Japan, the prevailing direction of propagation is 
somewhat westerly shifted from the due north direction having the 
shortest distance of the wave propagation, and the frequency distribution 
is roughly in inverse proportion to the distance from the Observatory to 
the land shelf. When the disturbance sources exist i i 
: Suis € sources exis 

B | ces exist in the Pacific Ocean, 
the prevailing direction is southerly, in whose direction the distance of 
propagation is longest, and the frequency has the tendency to decrease 


in proportion to the distance from the Observatory to the coast. This is 
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probably due to the position of the typhoon where the typhoon does not 
send swells to the south-easterly coast of the Kii Peninsula. The 
frequency distribution is found to be rather dense in the direction of 
Osaka Bay. These waves coming from that direction may be not 
considered to be generated at Osaka Bay separated almost from the 
open sea, They are supposed to be the waves propagated from Tosa 
Bay, because their periods are 
similar to that of the waves com- 
ing from the other directions. In 
the direction of Wakayama the fre- 
quency is minimum without any ex- 
ception in all distributions, This is 
the most interesting fact, and if we 
take up a pursuit of this fact, the 
origin of microseisms may be found 
soon unexpectedly. The frequency 


distribution shown in Fig. 18 was 


Wakayama Kii Peninsula 
obtained on the waves of the typical 


Rayleigh form selected among the Tosa Bay 
all waves observed by vector seismo- * covet ee 

2 Fig. 18. Frequency distribution of ar- 
graphs. This figure also shows that rival directions of microseismic waves 


all directions point toward the coasts having the pure Rayleigh type. 


with two exceptions and the frequency is minimum in the direction of 


Wakayama district. 


6. Progression of microseismic waves and 


their wave velocity 


For the purpose of investigating whether the waves of the Rayleigh 
type selected to make the frequency distributions are progressive or not, 
the observation by ordinary seismographs at the two points were carried 
out simultaneously with the observation by vector seismographs. The 
positions of the two points are shown by N and S in Fig. 11, and the 
instruments used for this purpose are NS-component seismgraphs with 
the same constants as those of the vector seismographs. The arrival di- 


rections were measured on the vector seismographs and the differences 
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Fig. 19. Differences of the arrival times of microseismic waves at the two points 
of observation versus the arrival directions. 


of the arrival times of micro- 


3 50 ° a seismic waves were read from 
e . the seismographs of the two 
40} a He points of observation. These 
ak oes Sines Se ae ag Ro ou relations are plotted in Fig. 
i o “ee i ne = ‘: ° 19. The three sinusoidal 

20450 100) «GCOS Sk curves in the figure show the 


Amplitude ratio, Aw/A, relations between the arrival 


Fig. 20. Wave velocity as a function of the directions and the 
amplitude ratio of horizontal to vertical 
component. 


differences 
of arrival times at the two 
points computed for the waves 
of the velocities of 3.6 km/sec, 3.0 km/sec and 2.4 km/sec. The points 
are scattered. But they are on the whole distributed between the two 
curves showing the velocities of 2.4 km/sec and 3.6 km/sec, while the 
velocities observed by tripartite measurements are widely scattered. The 
open circles in the figure show only the waves with the orbital motions 
of the typical Rayleigh form. These points are plotted along the curve 
of 3.0 km/sec. Fig. 20 shows the wave velocity as a function of the am- 
plitude ratio of horizontal to vertical component, which was deduced from 
the analysis of orbital motions. The wave velocity, especially for the 


waves of the typical Rayleigh form, is not scattered so much within the 


range of the ratio of 0.8 to 1.2. Accordingly the wave of the typical 


Rayleigh form may be considered to be the single wave, and there is little 


doubt that they progress with the velocity of about 3 km/sec. 
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7. Summary 


1) The microseismic waves are undoubtedly of the Rayleigh type. 
Even though vertical ground motions are small as compared with hori- 
zontal motions, it is in consequence of superposition of waves: coming 
from various directions. 

2) Vertical ground motions are nearly equal in mean amplitude to 
horizontal motions. The mean amplitude in NS component is slightly 
larger than in EW component, because the Rayleigh waves approach from 
the northerly or southerly directions in which the origins of generation 
of microseisms exist. 

3) The mean period of microseismic waves is nearly equal in each 
component, and it is about 4.1 sec. This period is approximately equal 
to half the mean period of swells in the westerly Pacific Ocean. 

4) There are the predominant period of 2 sec. beside of 4 sec, in the 
frequency distribution of the period. The 2-second microseisms are as- 
sociated with the local wind blowing at the time of the passage of cold 
fronts. 

5) The amplitude and the period of microseismic waves depend on the 
scale of the disturbance source and the distance of propagation of swells 
from the disturbance source to the coast near the station. 

6) The tripartite observation does not give the satisfactory results. 

7) The wave velocity is about 3.0 km/sec. 

8) The arrival directions observed by the vector seismgraphs indicate 
decidedly that the microseismic waves propagate from the coasts. The 
observation by vector seismographs is one of the most effective means 
to investigate the microseisms. 

9) All results of the writer’s investigations show that microseisms 
are generated by swells propagated from the disturbance source to the 


coasts near the observatory. 
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10. Free Oscillations of the Earth 


By Hitoshi TAKEUCHI and Masanori SAITO 
Geophysical Institute, Tokyo University 
(Comm. by C. TsuBol, M.J.A., Jan. 12, 1961) 


In view of recent progress in seismometrical instruments, it is 
of interest to calculate theoretically the periods of the earth’s free 
oscillations. This is a report of our study in this direction. 

Assuming that the density po and elastic constants 2 and pv of the 
earth are given functions of distance r from its center and denoting 
the (7,@,¢) components of displacement w and an additional potential 
P due to the earth’s deformation respectively by 


u,=U(r)S,(6, 9), y= V(r) eae 
up= (7) 9549) P= P(W)S,(6, 9), Se 


sin 0 a) 

where a common time factor e’’” is omitted and S,(0, ©) is a spherical 
harmonic of order », we can study the problem of the earth’s free 
spheroidal oscillations by solving three simultaneous differential equa- 
tions among U(r), V(r) and P(r). These three equations with homo- 
geneous boundary conditions to be satisfied at the free surface and at 
several discontinuity surfaces within the earth will determine the 
period T7=2z/c of its free oscillation as a function of the azimuthal 
wave number n. Making use of the earth model by Jeffreys and 
Bullen,” Alterman, Jarosch and Pekeris” calculated T for n=0, 2,3 
and 4. 

By using numerical integrations, we extended the calculation for 
the same model up to n=5, 6, 8, 10, 16, 24, 38 and 96. As is well known 
from previous studies,” the Jeffreys-Bullen model is not a very accurate 
approximation to the physical structure of the upper mantle. In view 
of this, we also calculated 7 for the Gutenberg model.” Since the 
main difference between the Jeffreys-Bullen and Gutenberg model of 
the earth is the existence of a low velocity zone around 150 Km 
depth in the latter, we may expect that a significant difference for 
the two models will be seen in 7 values for larger . Thus the 
calculation for the Gutenberg model is made only for n=16, 24, 38 and 
96. In making the calculations for »=16, 24, 38 and 96, we put a 
homogeneous crust of thickness 35 Km, density 2.8 g/cm’, compressional 
wave velocity 6.10 Km/sec and distortional wave velocity 3.535 Km/sec 
on both of the models, whereas we assume no crust in the calculations 
for n=0 to 16, For both of the models, the values of 7 for n=16 
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with and without the crust agree to the third significant figure, 
justifying the above way of handling the crust. The values of T 
thus calculated are shown in Table I. In Table I are shown also the 
values of 7 obtained by Benioff, Press and Smith® from spectrum 
analysis of records of the Chilean earthquake of May 22, 1960. 

In 1959, by using variational calculus, Takeuchi” calculated free 
periods of the earth’s torsional oscillations for the Jeffreys-Bullen 
model. Extending the calculation, we get 7 given in Table Il. In 
Table II, are shown also the observed 7 by Benioff and others.” 


Table I. Periods of the earth’s free spheroidal oscillations 


n | calculated T: min observed 7: min 
2 53.7 54.7 and 53.1 
oar" | 35.5 35.9 and 35.2 
at a 25.7 | 25.8 

5 19.8 | 19.8 
ony 15.9 | 16.0 
en 11.68 11.81 
100-5)! 9.68 | 9.66 

Jeffreys-Bullen | Gutenberg 

1620) 6.74 | 6.74 6.78 

24 5.05 5.06 5.10 

38 3.62 | 3.65 3.66 

96 1.658 | 1.679 | ee 


Table II. Periods of the earth’s free torsional oscillations 
ee 


n | calculated T: min observed 7: min 
2 43.4 42.3 
3 28.1 28.6 
4 21.5 28 
5 | LT oi: 17.9 
6 | 15.3 15.5 
fl 13.5 13.5 
8 Aa 12°38 
9 | 11.04 1G EAPAL 
10 10.17 10.38 


11 | 9.48 9.60 


From Tables I and II, we see that the calculated periods agree well 
with the observed values. Thus it is evident that the models used 
are good ones for representing the internal constitution of the earth. 

Furthermore, by comparing the calculated results with the obser- 
vations of mantle Rayleigh waves, we may see which of the Guten- 
berg’s model and the Jeffreys-Bullen’s is a better representation of 
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the earth. Calculating the phase velocity C of the mantle Rayleigh 
waves by means of the formula 


= 2na, 
(n+0.5) T’ oy 
where a is the radius of the earth, we get the relations between 7 
and C of the waves shown in Table III. 


Table III. Relations between 7’ and C for the two earth models 
eS ee eS eC 


Gutenberg Jeffreys-Bullen 
T: see C: Km/sec T: sec C: Km/see 
100.7 4.118 | 99.5 | 4.168 
218.9 4.749 | 217.1 | 4.789 
303.8 5.378 303.3 5.387 
404 6.00 / 404 | 6.00 


The group velocity U as a function of T as experimentally determined 

by Ewing and Press® is shown in Table IV. With assumed initial 

values C=4,00, 4.10 and 4.20 Km/sec for T=100 sec, the equation 
ae. CC=U) (3) 
dT Bie & 

is numerically integrated and the results are shown in Table IV. 


Table IV. Experimentally determined U and C 


T: see | U: Km/sec | C: Km/sec 
100 3.76 } 4.00 | 4.10 | 4.20 
150 3.66 | 4.159 | 4.824 | 4.493 
200 3.59 | 4.874 4.621 | 4,888 
250 | 3.58 | 4.681 4.984 | 5.872 
300 | 3.73 | 4,902 5.385 5.943 
350 | 4.00 | 5.185 5.766 
400 | 4.10 | 5.335 | 6.127 | 


In Table III, both Gutenberg and Jeffreys-Bullen models are seen to 
give C=6.00 and 5.38 Km/sec for T around 400 and 300 sec. Table IV 
gives C around 4.10 Km/sec as the initial value at 7=100 sec which 
is compatible with the above C’s around 7=300 and 400 sec. Thus 
in this respect the Gutenberg model, which gives C around 4.10 Km/seec 
at T 100 sec, appears to be better than the Jeffreys-Bullen’s. 
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Love Waves in case the Surface Layer is Variable in Thickness 


By 
Rydsuke Sato 


Geophysical Institute, Faculty of Science, The University of Tokyo 


1. Introduction 


Recently, the seismic model study has made 
remarkable progress and wave character- 
stics of various complicated structures have 
een studied experimentally by many authors. 
. Press (1957) investigated flexural waves, 
yhich are easily excited and propagated in 
lates having various shapes (thickness change, 
thology change, fault and scarp). A. TAKAGI 
1959) studied Rayleigh waves along the 
urface of a layered structure of variable 
aickness and reported many interesting pro- 
erties of the waves, for example, that when 
ayleigh waves are incident from the side 
f thicker surface layer to that of thinner 
arface layer, the transmitted Rayleigh waves 
haracteristic for the thinner layer can be 
asily observed even at points close to the 
‘ansitional zone, but that when Rayleigh 
raves are incident from the thinner layer 
de, the transmitted Rayleigh waves maintain 
haracteristics for the thinner layer even at 
oints distant from the transitional zone. 

It is important to treat such problems 
1eoretically. In this paper, Love waves 
ropagated in layered media of which the 
face layer is of variable thickness are 
vestigated by means of the method developed 


vo.1=A cos B2,wH-exp {—f1,w2—ik1,w("%—2%o)} , 
vo.2=A cos $2,w(z+H)-exp {—ik1,w(x—%0)} , 


here 


ra) / 
P2w=V Ri Bey 


SiN bi ee ’ 


1d «1,.y is a root of the equation 


tan B2,.wH=;5— 
Boy 


by D. S. Jones (1952) for some diffraction 
problems. 


§2. The case that Love waves are incident 
from the side of thicker surface layer 
to that of thinner surface layer 


We consider a layered structure with a 
surface layer of which the thickness is H for 
the range «>0 and h(=H-—6) for the range 
x<Q with the origin of coordinates taken 
along the interface between two media, with 
the z-axis vertically downward (Fig. 1). 


(9+4)-=H- 
Kiowel® 


Velocities of shear waves and rigidities are 
taken Vi, #1 in the substratum and V2(<Vi), 
f2o(<pu) in the surface layer, respectively. 

Let us express incident Love waves, omitting 
the time factor e~*’’, as follows; 


220, 


\ Recent 
0222 —H, 
Wo ki=o/Vi , k2=o0V2 , (2) 
rae, (3) 
2 


id if we write «1.v=/C:,v, Ci.w represents the phase velocity of Love waves of the N-th 
ode in a layered structure with a surface layer having the thickness H. 
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Denote the total displacements by 


V=V0,1+1 , 2 OF o>ey>—o, (42 
V=V0,2+02 , 0>z>-h, Co fe com (At 
V=V0,2+03 , —h>z>-H, e>0. (4 
Then boundary conditions are given by 
(i) 2=0, ~>"%>—00; 
vi=v pies (5% 
ESE: Ozer eat 
Ci) z=—h, 220; 
= Ov2 _ Os . 
V2=03 , aaeiea & ( 
(iii) z=—h, «<0; 
ous — Lay ye sin Be,wO-exp {—ik1,w(x—o)} ’ (5 
0z 0z 
(iv) z=—H, «£20; 
Ovs 
— =0) ; 5 
0z ( 
(v) #«=0, —h>z>—H; 
oO as = Ue i Ags cos B2,w(z+H)-exp {ir1,wao} . (5 
hy Chy 
v1 satisfies the equation of motion 
Oni Cos. Lewes 
aut! a2 Vi oF 4 
and v; (j=2, 3) the equation 
0’v; Ov; AE 070; (" 
ac? eee? Yor 


If we introduce a virtual resistance ¢(dv:/dt) («>0) to (6) and assume the time factor e-* 
equation (6) is reduced to 
0701 Orv1 
Ox" oz? 


+-3701=0 , (§ 


where k:’?=(*+iew)/Vi*, that is, ki has a positive imaginary part. Hereafter, we consid 
that _% ki>0 and ki=o/V: as the limiting case of e>0. vz and vs satisfy the equation simil 


to (8), substituting k2 in place of ki, and _* kz should be positive and [Ril <|Rel. 
If we write 


il oe , 1 0 7 
Vib, =e \ ‘ee Tm\_ etde=Vis(C, 2)+Vi(t,2),  CHe+in, ( 


then, according to a generalized Fourier’s theorem, if |v:|<A exp (n-#) as a—>-+00 a 
lui] <Bexp (y7+") as z>—oo, then Vi+(C, z) is analytic for 7>y_, Vi_(C, z) is analytic for n< 
and therefore Vi(¢, z) is analytic in the strip 7.> »>7-. In our case, since the displacemer 
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consist of reflected waves at «=0, transmitted waves and diffracted waves, we can take 
+= 7 ki and 7-=—F ki. 

Applying the Fourier transform in x to the partial differential equations for 1: and ve 
(cf. (8)), we find that 


a % 

aa US Z)—687 Vi, z)=0 , Bi wy C—k2, (10) 
da? o~ 2 ) oa aa aT ROR ESS 

aE V2(C, z2)—Be? V2(C, z)=0 , Bro=/ C2—R,? . (11) 


These equations have solutions 
ViiG, 2)+ Vi-(C, z) = Be-*:# , (12) 
Vos(C, 2) + Ve_(C, 2) =Ce-P2*+ De®a? , (13) 


where B, C and D are functions of € only. From the boundary condition (5a), we can write 


— Pat 7 Br Zeta T By 
(C= 2B: Be iD= yrs (14) 
Hence, 
V2i(C, z)+ Ve_(C, == (0 cosh f2z—7f: sinh f2z]B . (15) 


For brevity, we shall write Vii(C), V2s(C), etc. instead of Vis(C, —h), Vex(C, —h), etc..... 
' Differentiate (15) with respect to z and put z=—h. This gives 


Vo4"(€)+ Vo_"(C)=—[f2 sinh 62h +781 cosh B2h|B . (16) 


Eliminating B from (15) and (16), we obtain 


aed Bs cosh B2z—7 81 sinh f2z a, <% 7 
V2i(€, z)+ V2-(, 2) Bs FAO [Vew(C) + Vea"(O)] (17) 
F.(C€)=f:2 sinh B2h+7 8:1 cosh Beh , (18) 


or 
1 f: cosh P2h+7f: sinh Beh Vase! 
a a [Var'(6) + Ve-) 


(19) 


V2+(C)+ Vo-(C)= 


Next, we apply an operater Jal, te to the partial differential equation for vs. We get 
T Jo 


@ re eee Paine a 20 
ae Ve z)— Be? Var(C, Z) a Seal ax 1 Vig eee ’ ( ) 


where (0v3/0x)z-0 is given by the boundary condition (5e) but (vs)z-0 is still an unknown 
function. In order to eliminate the unknown (vs)s-0, changing the sign of € and adding the 
resulting equation to (20), we find 


[Veal 2)-+Vsu(—£, 2)]—Bs*[ Vax, 2)+ Vor(—6, 2)] = Fa C08 Bs. w2+H)-exp {irs,wate} , 
(21) 
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the solution of which is given, taking into account the boundary condition (5d), by PS 


ay! 21AK1,7 COS B2,.w(z+H) 
Vai(€, z) + Vai(—C, z)=E cosh B2(z+ 7) oe 2 8 


exp {irs _wito} . (22) 


As before, we differentiate (22) with respect to z and put z=—h. Then we have 


Var(O+ Ver —C)=E fe sinh 220-4 Sr Sin 62.99 on (ies wo} « (23) 


Been: 


Elimination of E between (22) and (23) gives 


21 An,wB2,w Sin B2,wd 
V 2a Cm x 


1 cosh £2x(z+H) 
Bz sinh $26 


2iAm,n Cos b2,w(zt+H) 
Von Oe 


Vike DEV KO = [ var@+ vers eo (ie. | 


exp {ir1,w2o} , (24) 


hence, for z=—h, 


2iAr:,wbe,n Sin B2,wd 
V/ 2 CH x 


V34(€)+ Vil -O= coth B2d [ vo } Va4'( C) exp (ies x | 


2iAni,w COS f2,wd 
/ on C—eaw 


exp {tK1 _wato} 7 (25) 


Now, from the boundary condition (5b), we can take 


Vo+(C)=Vae(0), = Va4"(C)=V4"(C) , (26) 


0 
and from the application of an operater Jz | dxe's* to (5c), we can get 
‘(a —coo 


iABs.w sin Bo, wd 


V2-"(Q)= WT ee 


exp {1k1, Xo} 4 YF (Cra) <0 5 (27) 


Using (19) and (25) to (27) and eliminating Vs,(¢) and Vs:/(€), we can obtain the following 
equation of modified Wiener-Hopf form; 


ile en ee : 1 re ; 
eee a [Vist(©) + Ve-"Ol + coth B:B-[Vee"(—0) + Va-(—O) 
+1 C)= Vil] SS eA xn Gea ne (28) 
where 
F,(€) =f» sinh B2H+ 78: cosh b2H (29) 


and F.(f) is given by (18). 
The factorization is made as follows; 


ee) 
sinh b26 F3(€) 


= hOKeOe (30) 
where 


K@=K(-D= 2 coat G1) 


Love Waves in case the Surface Layer is Variable in Thickness 23 


H.(€)=H_(—2) is given by (A6), L4(€)=L-(—£) is given by (A13) in Appendix and #:,n and 
Kam are Zeros of Fi(¢) and F.(f), respectively. K,.(0) is regular in _% (€+k:)>0 and K_(£) is 
regular in _% (€—k:)<0. Furthermore, we decompose coth 20/82 as 


a coth Bd= f+ £10), (32) 
where 
ft =f es (33) 
2k20 C+k2 n=1 pnd C+tp, 
and pf» is given by (A2) in Appendix. 
Insert (30) and (32) in (28) and rearrange to find 
K+) , tABs.~ K,(1,7) sin Baw - 
Eee Le 77 irene ey wre. MEAP CIT) 
i 1 V2 (Re) uch kot+ipn V2 (tpn) =P(0) (34) 


C+hke K4(k:2) n=1 Dn C+tpn Ki(ipn) 


in which V2_’(€) [V2’()= V2i(€)+ V2_"(O)] is, of course, a known function and is given by (27). 
Then, the left-hand side of (34) is regular in the domain 7% (€+:)>0 (C=m1,y is not an 
irregular point since the left-hand side is zero at this point) and the right-hand side is regular 
_ in the domain _% (€—h:)<0 and therefore a function, defined by the left- or right-hand side, 
is regular in the whole of the ¢€-plane, since the above-mentioned two half-planes overlap. 
Hence, if the left- and right-hand sides tend to zero as € tends to infinity in appropriate half- 
planes, both are zero from Louville’s theorem. If we assume (dv2/0z)—constant xu~\/? as 
x—>+0 on z=—h, then | V2+/(£)|<constant x |C|-'/.. From this edge condition, we can take 
that each side of equation (34) is identically zero. This edge condition, which is given by 
the physics of the problem, plays an important réle in the solution of this problem since it 
is concerned with the uniqueness of the result. 
Thus we can obtain 


Wwi@=— were t oa 1 Lo poke aoe (35) 


Ki(0) C—K1,¥ C+hke n=1 Dn C+ipn 
where 
walhe ae sin B2,wd-exp {ir1, ao} , (36) 
LV aiths) - I Ve ipa) 37 
aya Ge tee 2 Waki) © ete 


An infinite set of simultaneous linear algebraic equations for %, is obtained by setting 
a=? pa m=O, 1; 2,---) in (35): 


ECE ia ets ee, (38) 
K1,v—Re n=1 Pn 
[Ks (iPm)|?%m =i ko+tpm 14 > Ro+tpm Rattpny (39) 


k1,v—1Upm n=l pn DPm+pn 
- Furthermore, if we set C=—Az or —ifn, the following relation is easily found 


Ve!(—he) = — V2! (Re) , V2!(—ipn)=— V2'(ipn) . (40) 
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We use the Fourier inversion formula to find 


masge| Val ae*edl, K£%kh>c>—-Fh. (41) 
/ on 


ie—oo 


Therefore, from (17), we have 


1 (i+ 1 B. cosh foz—yfi sinh Pez yy, Ce-i 
ae eee e s*de . 42 
ie Ti. Bs FD) segs 4 
Next, vs in #>0 is given by 
—1_[""y.@, de-eedt (43) 
m= sec\ ail, Ze , 


but, since Vs,(—C, z) is regular in the lower half-plane, 


[Vel zje*s*dl=0, «>0, 


to— 


hence, we can also write as 


v= Te. Val, 2)+Vaul—6, aleedl, —-w>0. Soin 


From (24) and (27), we find 


ice —oco 


Vat D4 Vatet, Jo oe Dy, 


()+ V2!(—O)] 


Be sinh B26 
iA x 1 1 t 
tie Oe Benle+H) —— st et EXP {tK1,wao} . (45) 
Then vs in x>0 is given by 
Vs=Us+ A cos B2,w(z+H)-exp {iri,w(a+a0)} , x>0, (46) 

- . b (#9) cosh Siz+H) 
e————— y Vad irae eer 
: Til... are mmc iS (47) 


First, we will investigate the transmitted Love waves, i.e. ve in w<0. In the region 
0>z>—h and « is large negative, we can close the contour in the upper half-plane. V2/(C) 
has only a pole at €=*,y in the upper half-plane. From (35), 


» tABe.w . x : 
[V2"(C)(C—m1 gm, y= — Aree Sin S2,wd-exp {ik1,wXo} , (48) 
then, using the relation (3), it is found that the contribution vz. from this pole gives 


v2i1=A sin Bz, yd- B2,.w COS B2,.wz—7B1,w Sin B2,wz 


Baw sin B2,wh—7YBi.n cos Bsaiehe exp) {—1K1,w(a@—20)} 


=—A cos 62, w(z+H)-exp {—ik:,w(a—a0)} , 20, (49) 


which exactly cancels the incident Love wave v2. 
Denoting the zeros of F2(¢) in the upper half-plane by 2m (m=1, 2,::+), these poles give 
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v2=—/dxi > ei B2,m COS Be,m2—7B1,m Sin Be.mz 


V2! (Ke >) “exp {—ike mic} 


Ho Bo,m [dF 2(C)/df]e- Kom 
ae = Kam C2,m cos Bs, m(z-+h 5 
V2xt 2») Bia Ges =) Sa zon, (k2,m)*exp {—ike, ma} , £<0, (50) 


where 
Bi,m=/k, 2—Re » Bam=vV kee, 2) 


C2,m=@/k2,m: phase-velocity of Love waves of the m-th mode propagated in 


the layered structure with a surface layer having thickness h, (91) 
Uz.m: group-velocity of Love waves corresponding to the phase-velocity 
(Ons 


These waves represent the transmitted Love waves. 

Next, we will evaluate the reflected Love waves, i.e. vs in x>0 (v2 for large x is identical 
to vs). For the large positive x, we can close the contour in the lower half-plane. In (47), 
zeros of $2 sinh /20 and sinh §2d are given by C=—kz and C=—ipm (m=1, 2,---). At these points, 
however, V2’(¢)+ V2/(—C)=0 as shown in (40). Therefore, these points are not poles. In 
the upper half-plane, V2’(¢) has poles at C=—K1,m (m=1,2,---) and V2(—f) has a pole 
€=—m,y. At C=—nr1,y, 


[Ve(—O\C Hes ees = = 14h sin Br, w-exp {irra} , (52) 
then, the contribution v3,1 from this pole is 
Vsi=—A cos Be,w(z+H)-exp {ix1, w(a+ xo} -) x>0 5 (53) 


which cancels the second term in the right-hand side of (46). The contribution v3,2 from 
the poles €=—v*«1,m is easily obtained as 


Fo.2=a/ Tei Se COS Be mH) yy (004 1 m)lga-e, EXD fier mt}, > 0, (54) 
m=1 [2,m SiN B2,m0 


which represents the reflected Love waves. 
§3. The case that Love waves are incident from the side of thinner surface layer to that 


of thicker surface layer 
Write the incident Love waves as 


Vo,1=A COs Ba arh-exp {(—Br.az-+ ies, + a0) : z20, \ ee (55) 
vo,2=A cos B2,w(z+h)-exp {ire w(a%+%0)} , 02>2>-h, 
where 
Biw=V 5, p—Re> — Ba.w=V Rey s > (56) 


and «2,v is a root of the equation 


“a 


tan Be,wh=p-s . (57) 
Bz 


es 


Displacements are taken as 
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| v=vVoitMN , 220, o> 17> —oo, (58a) 
V=Vo0,2t+0V2 , 0>z>-h, Copii — C8) 5 (58b) 
V=0s ; —h>z>—-H, “20, (58c) 


then, the given boundary conditions are 
(i) 2=0) arr —o; 


-_ Ovi _ Ove (59a) | 
ee i az or 
(ii) z=—h, 220; at $e | 
: v Dv 
Vs—V2= A-exp {ike w(a+Xo)} - ee - (59b) 
(iii) z=—h, w<0; 
Ba eaiys (59c) 
0z 
(iv) z=—H, «20; 
Ge (59d) 
dz 
(v) #=0, —hoz2—H; 
a (59e) 
hy 


Proceed in the same way as before up to equation (19). For V3:(€, z), since (0v3/0%)z-.=0 
in this case, we have, instead of (24), 


Vasll, 2) 4 Var(—, z= Shh Dy, OL VED). (60) 
Bs sinh B26 
Hence, at z=—h, it is reduced to 
Vs4(C)+ Vsi(—0)= ; coth B26[Vsi’(€)+ Vsi’(—f)] . (61) 


Use the boundary conditions (59a, c) to find 


Vas"(C)=Var'(S) , Va+(0)—Vas(6)= ae oe exp {tke wo}, F% (E+kew)>0, 


(62) 
Ve2(e\= 0" 
Thus, from (19) and (61), we obtain the equation, corresponding to (28), 
1 eden Fe 1 
oe a Var") + coth B20 Vas"(—0) +[Va©)—Var(—0)] 
iA 1 A 
on Ea exp {ix2, ao} =0 k (63) 


According to the same factorization as before, it is found that 


K,(6) {Ad’. bee eee 
ee Vay (C) + 
C+he HG) + V/ or LAO hee exp {7K2,v Xo} 
1 V4."(Re) kotipn Va+"(ipn) 


ae 


i 
C+hke K4(k2) n=1 Pn C+ipn Kx(ipn) =P_(C) ’ (64) 
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and this corresponds to the equation (34). Thus, applying Louville’s theorem, we have 


& ke 1 1 See. ip 
Vs "O=— wet eae MiB See 
: K4(6) abr Clb n=1 pn €+1pn ay (oy 
where, in this case, 
a Ao ko+k2,w ex {i \ 
On K4(2,7) ae a ae 66 
pe tleg Vs4/(Re) 7 _ 1 Vs4"(ipn) oo 
W K(k) ’ “ W Kx(ibn) 
Besides, the relations similar to (40) are also obtained as 
V34/(—k2) = — Vas" (ke) , V34!(—tpn)=— V+" (ibn) . (67) 
Now, from (17) and (62), we can write 
i Be cosh Boz—7 1 sinh Boz 
V2(C, z)= i(O- 68 
2(C, 2) Bs FAO) 3+/(C) (68) 
hence, the displacement v2 in 02z=—hA is given by the following integral; 
__1_ (*t" 1 Be cosh fez—rfi sinh 22 1 yey iene 69 
a Tiel... , as Vorl(O)-e-idl . (69) 
v3 in x>0 is also given by (44). Using the relation (60), we have 
1 (*t- 1 cosh f2(z+H) HO) Vas! ~igeg 
a = i d So). 70 
= vale Be sinh f26 ys Or - ( Oe : : 


The transmitted Love waves can be obtained from the residue of integral in (70). For large 
positive x, we can complete the contour in the lower half-plane. C=—k, and C=—tpm 
(m=1,2,---) are not poles of the integrand according to the relations given by (67). Poles 
inthe lower half-plane are only —«1,m(m=1,2,---), which are zeros of K.(¢). Thus, the 
transmitted Love waves are given by 


ne oa > 608 Bam ED TY, (OC + 1 m ees "EXD (er ne OF 


m=1 B2,m Sin Be, (71) 


Bem=V hte, A 


which is of the same form as (54). Poles of integrand in (69) are given by «2, (n=l, 2,-°*); 
which are zeros of F:(¢), in the upper half-plane. Thus, closing the contour in the upper 


half-plane, we finally obtain 


1 §B2,m Cos B2,mZ—7P1,m sin B2,mz Vou e2,m) -exp {—2K2, ma} 


ey ca [AF (C/E lees m 
ek oan (Ou Cos Bz,m(z+h) , = 
ie 2m Se ee V oul (hain sCXD {— 182, me} x<0, (72) 
V2nt 2 B2,m ae ) B2,mh Bs 


as reflected Love waves. 
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§4. Approximate solutions for the case of small /2d 

In order to evaluate the reflected and the transmitted Love waves numerically, we must 
solve the infinite set of simultaneous linear equations, for example as (38) and (39), and even 
if we could obtain the solutions of the above linear equations, we must also calculate K+(#2,m) 
in V2!(r2,m) (see (50)). Besides, to obtain the exact values of K+(xz,m), we must get all zeros 
of Fi(€) and F.(¢), defined by (29) and (18), for all frequencies o. 

Here, we assume that k2d=k2(H—h) is very small, i.e., the difference in thickness is very 
small compared with the wave-length of waves under consideration. 

(a) The case when Love waves are incident from the side of thicker surface layer to that » 
of thinner surface layer (Case I). 


When € is real, H:(¢), defined by (A3), can be written as 


- 1/ = 
Hy Q=[SRS" | exel-2O +4 ZCoe—on)], Ga=tan-' Cbs). (73) 
2 n= 
Hence, if we omit the term of O(%20) compared with 1, we find 
H(k2)~1, (74) 
then 
K.(k2)~1, (75) 


to the approximation of the same order. 
Next, when C=ipm~1/6m=1(mz/6), since 


ip (1 eegpe ee ) 7, be <1n(1—C—In m) (76) 


T T 


(see A7)), we have 


ELA) = 1 ( DnOn+ pmOm) exp {— pmOn—y(ipm)} ~ Il (1+ ee (1—C—I1n m) 
n= n=1 n 


m”e-™ 
Sas OT (77) 


The remaining terms of K,(ipm) are 1+O(k26). Thus we can write 


dw le 


K(tpm) “Ta +m) . 


(78) 


When m is very large, K,(ipm)~1/\/2zm. Omit the term of O(f26) in (38), (39) to find 


[1+ {K(he)}*]y0~2——2ie__, 
ko—k1," 

Ey. : = mM (79) 

{K+(tpm)}°y¥m~—i—y— > an 
n=1 M-+N 
Hence, from the above approximations, 
Ba ike 
ae Ro—k1 N 

m2™e-2m Kiw © (80) 


es a ES m 
(ray mee a 
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In order to obtain V2/(f) given by (35), we must know the value of Kyi(k1,~). We can 
easily get Ky(«1,v)~1 from the same approximation as K,(kz), taking into account that 
k1,wO =(K1,w/Rz) (Rod) <(Ri/k2)k2d<k2d. K4(x2,m) nearly equals to I (wore). 

Thus, V2’(«2,m) in the right-hand side of equation (50) is given by 


i 7 


Ke,m—Ki,w  — Rotke,m n=1 


, ABs.w Rethem . 2 R : 
V2 Cham) oe ee sin B2,vd exp Ges. waa} | 


(81) 


which is of the order (k26) for m#N. When m=N, i.e., the modes of incident and trans- 
mitted Love waves are equal, we can obtain the solution as follows: 
Differentiate the equation 


SLE es 
tan EO hey ct ee A 82 
AS a f V kot —n? (82) 


with respect to J, keeping ki and k2 constant. This gives 


Whe — bal (E 
ante =t) , (83) 


where C=o/x and U=do/dx. Hence, for the same mode and the same frequency, 


fi, y eet (Eb =) One (84) 


Inserting this relation to the expression of V2’(2,w), we find 


oe ee aE BS 8 Be To a oe lO). (85) 
/ 2x ken _- =i) ti) 
Tae 


Thus, (50) can be finally written as 
v2,2~A cos Be, w(z+h) exp {—ixe,we+ ins, wato} [1+ O(F29)] . (86) 


This is the transmitted Love waves for small (k:6). It must be noted that all the waves 

with different modes from the incident Love waves are of the order (20). 
Next, we will estimate the reflected Love waves to the same approximation. From (A17), 
K(0) Laks) Il’ (C+41,n) (87) 


(C+K,,m) ~ Ay(6)(C+K2,m) n=1 (C€+K2,n) ‘ 


where II’ means an infinite product except m=m. According to the same approximate 
procedure, we have 


[ K.(6) | oe BAe ee (88) 
(C+41,m) $=-Kk1,m %2,m—K1,m Ki,m ae ) 3 


[Vo2"(O(C-+41,m)]¢-—«; m in the right-hand side of (54) is found to be 


A Ki,m Ro? +k1,K1,m ce 


ee 1) 6 sin B2,wd-exp {ik1, wo} . 
[V2"(C)(C+41,m)]¢=—«1m ee pie at PS ee ie ) sin 2,7 p 1, Vo 


(89) 
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Hence, we can obtain, as the reflected Love waves, 


032A Ay CoS B2,m(z+H) x ad exp {2K1,me+2K1, WX} . 


Bo,wh2,mEl  «1,09-+h1,m 


fi) ko? +k1,wK1,m (a2 =i) sin Bowd vO 
U1,m sin f2,m6 
(90) 


When m=N, it is reduced to 


~ ra) 2 2 
V3,2~A Cos B2,w(z+H)- ke ee N 


asi! Hey, w(% + %o)} « 91 
2H ha?—«, 3\Uw )exe parasea oa 


(b) The case when Love waves are incident from the side of thinner surface layer to that 
of thicker surface layer (Case II). 

We can evaluate the waves for this case by the same procedure. From (65), the equations 
corresponding to (80) are 


eee 
be: Retkew ’ (92) 
m2™e-2m ik2. oo m 
‘Theme bale ie eee 
Hence, 
GAO? Bot —=miakew kta Cal am , 
[Ve +41, gees mg aL 7 ( a 1) exp {ire wa} , (93) 
which is of the order (k:2d) only when m=N. For m=N, (93) can be reduced to 
[Vas'(OC-+ ew gant. we— EO By 2 exp {ire, woo} . (94) 
Then, we have, from (71), 
vs~A cos B2,w(z+H)-exp {ik1, we +ike, wXo} 5 (95) 
All the waves of other modes are of the order (£29). 
To the same approximation, 
LAO ke*+-k2,wKke,m ; 
Vou! sin) ee ALE 
3+/(K2,m) est eT exp {iK2, Ho} (96) 
Putting this into (72), we finally obtain 
mo ke? tke, NK G ' z 
n~—A Ss m H = 2 Lt) BY us == 
Sco: Bs. m(z+H)— Wer ene 1) exp {—ire,ma+ine, woo} . (97) 
When m=N, 
Pa: 4 6 ke +k, Gs. : 
vase — Acos Be. wl2th): 5 pe “(Gs 1) exp {—its,x(a—m0)} . (98) 


§5. Ratio of energy flux of transmitted or reflected Love waves to that of incident Love 
waves 


We could obtain the transmitted and reflected waves for small (k2d). These waves are 
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summarized as follows: 
(a) Case I. 
Incident waves; 


v0.1=A cos §2,.wH-exp {—61,wz} cos [wt+m1,7(2—20)] , z>0, \ So 99 
5 ml). 
vo,2=A Cos B2,w(z+H )-cos [wt+e1,¥(2—ao)] , 02z22—H, a 
Transmitted waves; 
Vt,1=A cos Bo.wh-exp {Bs wet cos [wt+-r2,we%—K1,vXo] , 220, \ <0 (99b) 
Vt,2=A cos B2,w(z+h)-cos [ot +2,~e—K1, wo] , 02z2-—h, 
Reflected waves; 
Uri=A Di am cos B2,mH-exp {—f1,mz} cos [ot—Ki,me—ki,wao], 220, 
‘ WW) (99c) 
joa Al = am COS B2,m(z+H)-cos [wt—K1,me—K1,w2X0] , 0>z>—-H, 
ei eO Ree + 1, wK1,m Cin sin Bo wo 
—_——-= 1)= eared ae 
Be,wh2.mll ki wekKi.m Ge Jee B2,m0 ae 
(b) Case II. 
Incident waves; 
Vo,1= A’ cos Ba.wh exp {—B w2} cos [wt—r2,w(a+20)] , 220, \ o<0. (101a) 
Vo,2=A’ cos f2,w(z+h)-cos [wt—K2,v("%+20)] , 02z>-h, 
Transmitted waves; 
vt,1= A’ cos Bs. H exp {—B1,wz} cos [ot—ei,we—re,yto], 220, \ eo: (101b) 
Vt,2= A’ cos 82,w(z+A)- cos [wt—m1,wu—ke2,w%o] , 0O>z2>—-H, 
- Reflected waves; 
Ur 1A! Dy Am! COS B2,mh-exp {—B1,mz} cos [wt-+K2,me%—K2,weo], z>0, 
m=1 B08 (101c) 
Vr.2=—A’ > am’ cos B2,m(z+h)-cos [ot +K2,ma—K2,wito] , 0>z>—h, 
m=1 
eee sae 1) ; (102) 
Bo alt Ke,nt+Kk2,m U2,m 


Since waves in cases I and II are of the same form as seen from the above equations, we 


will proceed as follows only for waves in case I. 
The rate of energy, which passes the plane perpendicular to v-axis, per unit time and 
per unit area, is given by 


FS p--—, (103) 


— being taken when waves are propagated to the direction of positive « and + to the 
opposite direction. Hence, we obtain 
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Fo.1= p10 Ates,w Cos? B2,wH- exp {—261,wz} sin? [ot-+e1.1(a—a)], 220, \ a>0. (104a) 
Fo.2= po A’e1,w Cos? Bo. w(z+H)-sin® [ot-+«1,x(2—a)] , 0>z>-H, 
Foe pate, coe Bz,wh-exp {20 2h sin? [ot+-«2,w%—K1,w%o] , 220, \ a<0. 
Fi,2=p20 A®k2,w COS? B2,w(z+h)-sin® [wt+K2,7%—K1,wHo] , 0>2>—h, (104b) 
F,.1= 0A? >) Am?k1,m COS? Bo.mH-exp {—261,mz}-sin? [ot—K1,m%—k1,nXo] , 220, ba>o 3 


F,2= p20 A? >) Gm?k1,m COS? B2,m(2+ H)-sin® [ot—K1,m%—K1,wH0] , 0222-4, 
(104c) 


jens x 
The energy flow for a certain period T=2z/o is, from F =| F dt, found to be 


0 


Fos= mp Ate w cos? Bs.wH-exp {—26:,w2} , \ ae (105a) 
Fo.2=7p2A?K1,~ Cos? B2,w(z+A) , 
Fia= np APra,w Cos? Pa.wh-exp {—2B1,v2} , \ 2<0. (105b) 
Fy,.2=7 2 A?K2, w CoS* B2,w(z+h) , 
Frist A? > Gn kins cos? Bs.mH-exp {—2B mz} s \ x>0 k (105c) 
Fy, .2=Tp2A? Di Am*k1,m COS? B2.m(z+H) : 
Thus, the total flux of energy across the vertical plane of unit breadth is 
na o = i Ki,w tan B2,.~H 8 ~ sec? Bow 
"ToP jing | ol'ot o. Fo,2 dz=—mp2 A*® = = aay wH—-, 
[ ] \ 0,1 zt | 0,2 dz 9 Tle Doe cca ht [1+ gat Se fot 
As Arfex Baw for incident Love waves, __ (106a) 
2 Ki,w = 1) 
Ui, 
wal » baw Be.wh : 
[T.F .Jirans. =—7p2 A for transmitted Love waves, (106b) 
2 Ko, ay 
U2,N 
and 
rTP leas, a Tyu2A* >) Am Bam BemH for reflected Love waves, (106c) 
2 els CS =) 
U1,m 
Hence, 
: Ge Pu (oe -1) (G2 1 
[TF Jerans. 1. Bax h Cie ee Gi N V2.2 Ui. (107a) 
[T.F .Jinc,  ke,7 Breas (fx =} H Civ (ene =) ie =) i 
U2.v 2" U2." 
and 
¥ G ng =i Cre se} Cin 1 
[T.F.]reri Sant Ki.w ee U1 9 = x os Ve Ui.w 
(al) Beiine Ki,m B,2 le m =) (Chie (3 mee! ) E m =) 
U; ™m Vv Ui m 
(107b) 


When m=N, i.e., the mode of reflected Love wave is equal to that of incident wave, it is 


Love Waves in case the Surface Layer is Variable in Thickness are 


easily found 


Rea Sa ee = 
Pick. lina, 


We. 


§6. Numerical examples 


(108) 


Fig. 2 shows phase- and group-velocities of Love waves (lst mode) in the layered 
Structure of #4i/4e=2 and Vi/V2=4/3, and [T.F.Jine./37p2A* for k2H is shown in Fig. 3. We 
must remember the assumption that kd =o(H—h)/V2 is very small. If we take 6/H=a, then 
k:H=k:6/a and k:h=k:H-(h/H). Hence, when a is not small, we can discuss only the waves 
for very small kxH. For an example, if h/H—0.8 or 6/H=0.2 and k:6<0.01, only the waves 
in the range kzH<0.05 or kzh<0.04 can be discussed. The smaller the value of h/H, the 


smaller the range that we can discuss for. 


And for the small values of 6/H or for the 


small values of keH, the value of [T.F.]Juans./[T.F.Jinc. is, of course, nearly 1. 
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Fig. 2. Phase- and group-velocities of Love waves 
under considerations. 
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Fig. 4 shows {[T.F.Jren./{[T.F.]inc,} X {1— 


(h/H)}~* for the reflected Love waves of the 
first mode, i.e., a?x{1—(h/H)}*. From this 
figure, we can estimate [T.F. ]irans./[T.F.]Jine. to 
some extent. 


§7. Conclusion 


We studied the problem of propagation of 
Love waves in the layered media with a 
surface layer of variable thickness. It is 
difficult to see the exact behaviours of Love 
waves in such media, because we must solve 
an infinite set of simultaneous linear equations 
and calculate phase- and group-velocities of 
Love waves for infinite numbers of modes, 
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as stated at the beginning of § 4. 
We derived approximate solutions for the case that the thickness difference of surface 


layer is small compared with the wave-length under consideration, leaving the solutions for 
larger thickness difference in the future. 

If the incident Love wave is of the N-th mode, only the transmitted Love wave of the 
N-th mode predominates and the waves of other modes are small. But the amplitudes of 
reflected Love waves are of the order (thickness-difference/wave-length) for all modes. 

The mean group-velocity Um for the total epicentral distance A is, of course, given by 


A/Um=A:1/U1+A2/U2 , 


for longer periods, as illustrated in Fig. 5. The period for Airy phase, of course, differs 
from that for a uniformly layered structure. 

The method used in this paper can be applied to solving the problem of propagation of 
Love waves in layered structures in which the boundary has a step form, as shown in Fig. 
6, with some modifications. 
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REEENDIXGa Decompositions at Sune ; ROT OK-© [ef. (30)] 


According to the infinite product theorem (E.C. Tircumarsu (1939), p. 113), we can 
write 


sinh B26 _ x 
Ed dh batt + Cant =H) , (Al) 
where 
Pndn=(1—h2?bn?)/? = —i( Ron? —1)/2, 
On=0/nn A i oo 
We decompose H() as follows; 
A (C)=H(C)H-(£) , H1(0)= Il { Pudn¥ 128n} exp {+720dn Fy} , (A3) 


taking x(C) as an arbitrary function. When € is very large 
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exp {ic — Ke) | 


H.(C)~ Il (1—i€dy} exp {tCdn—7(0)} = —_ z —, 
i) 


where C is the Euler’s constant and /(z) the Gamma function. From Stirling’s formula, 


(A4) 


we have 
$0 nf OER Ts, 0 @)) CO 
ra i =) / 35 exp {- qitee( 1 in ea (AS) 
hence, 
ey esl x. C6 ‘Cede er 
HC) S, 266 exp es i exp { th = 1—C—In =) +3 20 } * (A6) 
Thus, if we choose 
Des) Co\ 66 
(C0) ma C in) (A7) 
then, |7.(C)|~|C|-"” as |fl—-0o. Since 
AC) aCe (A8) 


H_(£) has the same behaviour for |C|—°co as A.(C). 
If we write all zeros of Fi(O) and F.C), given by (29) and (18), as tein and +ke.n 
(n=1, 2,3, ---), we can take 


PC) Se (Sei) GO) 
SWI 
Fae n= -1(C?—r, 2) #) Gl (Cj oo 


where 


G=FO/MC—-), GO=HO/MC—,), (A10) 


and both functions have no zero. Furthermore, we decompose L(¢)=G:(C)/G2(C) as follows; 


Gi(C) a ee 
L(O=— = Li (C)L-(O) , All 
(C) Gall) C)L-(C) (All) 
then, L.(©) is given by the integral 
a Li@=s5 | In LW) ay (A12) 
2 r w—C 


where we take the contour /’ as shown in Fig. Al, according to the definition of the upper 
sheet of the Riemann-plane (AF f:>0,.A f2>0). Hence, 


In Fw) tw a i In F(w) dw 


In L4(Q)=; : : | 
201 Jr 


w—C 2Qni w—C 
=e ee du A ee (A13) 
t So u—10 ZN WaRG 
where 
t/ wth? Cos V utt+heP tv ki —u? cos V ko? —w®H (A14) 


a, ; tan gy= 
tan di= V et hesin V uke J ke—u2 sind k2—wH 
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JOG a! 


and 2, %2 are given by the equations similar to (A114), taking h in place of H. And we 


have 
LA = Lito). (A15) 
From {(A1) and (A9), if we write 
__ PFs) KiG= Li(Q) LO) ll (C=ga ; (A16) 
sinh $25 F2(C) eo caer H(C)H-(C) a=1 (C—e, & 
we can finally obtain 
KO SK G2) 1 ee (A17) 


and |K+(¢)|, |A-()|~I¢1'? as |g|—00. 
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1. Introduction 


In the previous paper (Okano 1961) it was shown that microseismic 
waves are generated by the swells propagated from a disturbance source 
to the coasts near an observation point. Especially their arrival directions 
deduced from the orbital motions of ground which were recorded by the 
vector seismographs, made this inference reserve the writer's entire con- 
fidence. The further investigation was made by carrying out the obser- 
vations at the Aso and the Yura stations. The former station is located at the 
nearly central part of Kyushu where is surrounded by seas on all sides, 
and the latter is located at the coast of Wakayama Pref. in the southward 
direction at the distance of about 100 km from the Abuyama Seismological 
Observatory. Both stations appear suitable to investigate furthermore the 
precise location of generation of microseisms. The frequency of arrival 
directions of microseismic waves obtained at the Aso station is very par- 
tially distributed and these partially distributed directions seem to be 
closely associated with the slope of continental shelves. The frequency 
distributions obtained at the Abuyama and the Yura stations also support 
this close relationship. Furthermore the writer intended to investigate 
whether microseismic waves are generated at the regions close by coast 
lines or at some distance from coasts. The amplitude of waves coming 
from different directions supplied some information regarding that pro- 
blem. But the sufficient elucidation of this subject must be expected by 
further investigations. 

When the typhoon no. 14 in 1960 passed through on the south ocean 
off Japan, there were appreciable differences of the periods among the 
microseismic waves coming from different directions. This result is ex- 


plainable, if the periods of swells increase with an increase of distance of 
propagation. 


2. Arrival directions of microseismic waves at Kyushu 


With the same method used at the Abuyama Seismological Observa- 
tory, the arrival directions of microseismic waves were deduced by ob- 
serving the orbital motions of eround at Kyushu. The station where the 
observation was carried out is the Aso Volcanological Laboratory, located 
nearly in the central’ part of Kyushu. The instruments used are the vec- 
tor seismographs adjusted to the following constants ; T:=4.0sec., T2=6.0 
sec., hi, he=near critical and Vmax.=about 2800. Some of the orbital mo- 


tions recorded are shown in Fig. 1. It is found that they are appreciably 
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Fig. 1 Seismograms recorded at Aso station. 


disturbed and complicated in comparison with those recorded at the Abu- 
yama station, because of superposing by various types of volcanic micro- 
tremors (Sassa 1935). Notwithstanding the fact the waves of typical Rayleigh 
type, disregarding the slight deformation, are found so frequently that the 
frequency distribution of arrival directions can be obtained without dif- 


ficulty. Fig. 2 shows the three distributions obtained from the microseismic 


storms generated by the typhoon no. 6, no. 8 and no. 18 in 1960, and 


Fig. 3 shows the respective weather charts at the times when those dis- 
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tributions were obtained. When 
those distributions were obtained, 
the typhoons regarded as the dis- 
turbance sources were on the south- 
westerly ocean far off Kyushu. Most 
of the arrival directions do not point 
toward the centers of the typhoons, 
though there are a few pointed to- 
ward the south-westerly coast where 
considerable parts of swells arised 
in the typhoon areas are supposed 
to be sent, judging from the location 
of the typhoons. A great part of the 
frequencies is concentrated in the 
south-easterly direction (the Hyuga 
Sea), and a small part of them indi- 
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Fig. 2 Frequency distributions of arrival directions of microseismic waves 
observed at Aso station, 
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Fig. 3 Weather charts at the respective 
times when Fig. 2 was obtained. 


cates the north-westerly direction 
(the Genkai Sea), though the mi- 
croseismic storm on July 25 in 
1960 when the cold front passed 
over the Genkai Sea distributes 
considerable frequencies in that 
direction. And a little rests of 
the frequencies are distributed 
only in the direction of the Satsu- 
ma Peninsula. Previously the wri- 
ter deduced from the arrival direc- 
tions of microseismic waves at the 
Abuyama station that there are 


some districts where microseismic 


waves are frequently generated and some other districts not frequently. 


These observations at Kyushu supported more firmly this inference. 


Vertical components, as shown in Fig. 1, are small in amplitude as 


compared with horizontal, whereas the seismographs at the Abuyama 


station recorded the microseismic waves with similar amplitudes in two 


components. 
the Aso station is founded. 


This may be due to the property of the ground on which 


3. Arrival directions of microseismic waves 
at the Wakayama Pref. 


‘To obtain further information on the generating district of micro- 
seismic waves, the arrival directions were observed at the Yura observing 
station of the Wakayama Pref. by the vector seismographs used for the 
same purpose at the Abuyama station. 

When the seasonal wind stirred up the Japan Sea on Dec. 6 in 1960, 
the observation was carried out for about five hours. The seismographs 
recorded frequently the orbits of the considerably typical Rayleigh mo- 
tion, having the amplitudes in vertical component similar to those in 
horizontal, as recorded at the Abuyama station. An example of those 


orbital motions is shown in Fig. 4. Fig. 5 shows the frequency distribution 
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Fig. 4 Seismogram recorded at Yura station. 


of arrival directions read from those 
orbits. Notwithstanding that the 
Yura station is far more distant 
from the coast of the Japan Sea 
than that of the Pacific, almost all 
arrival directions point toward the 
Japan Sea. And the pattern of the 
distribution is very similar to that 
of the Abuyama station, that is, the 
microseismic waves do not come 
from the easterly and westerly direc- 


tions in which there are no coasts 


of the Japan Sea close to the station. : ta ; 
P Fig. 5 Frequency distribution of arrival 
The direction indicating the maxi- directions of microseismic waves ob- 


mum frequency does not coincide served at Yura station. 


with the northward direction in which the station is nearest to the 


Japan Sea, but deviates somewhat westwards. We may recognize the 
difference of about 10 degrees between the deviation of arrival directions 
at the Abuyama station and that at the Yura station. The existence of the 
difference appears to support the writer’s opinion on the origin of 


microseisms. 


4. The microseismic storm observed at the 
Abuyama station on Aug. 20 1960 


The strong microseismic storm was raised when the typhoon no. 14 
passed over the southerly ocean of Japan. The vector seismorgraphs 
recorded the orbital motions of ground in horizontal and only one 
vertical (UD-NS) planes, and simultaneously an ordinary seismograph 
with the same constants as the vector seismographs recorded the ground 
vibration in NS component. The arrival directions were read from the 
former and the wave period from the latter, because the vector seismo- 
graphs are hard to give the wave period so exactly as read from an 
ordinary seismograph. And if we intend to select only the waves of the 
typical Rayleigh form, the orbital motions in the horizontal and only 
one vertical (UD-NS) planes may supply the sufficient information for 


the arrival direction. Of course, if the microseismic waves come from 


p 
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Fig. 6 Frequency distribution of arrival directions of microseismic waves 
observed at Abuyama station on Aug. 20 1960 and location of the typhoon. 
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the easterly or westerly direction, the orbits in another vertical (UD-EW) 
plane are necessary to detect such waves. But the writer could save the 
necessity of those orbits, because he did scarcely find the east-westerly 
linear orbits with comparable amplitudes in the horizontal plane. The 
frequency distribution of arrival directions is shown in Fig. 6. The arrival 
directions do not point toward the center of the typhoon, and the direc- 
tion of maximum frequency deviates considerably eastwards from the 
south, differing from that obtained previously. This may be ascribed to 
that the typhoons were westwards from the Kii Peninsula in the previous 
case and eastwards in this case and hence the swells propagated from 


the respective typhoon areas were interruped by the Peninsula. 


5. Consideration on the district where microseismic 


waves are generated 


The frequency distribution of arrival directions of microseismic wave 
having the typical Rayleigh motions is shown in Fig. 7, including the 
new data obtained by the observation at the Abuyama station. Examining 
these distributions obtained at the 


Abuyama, the Aso and the Yura 


Noto Peninsula 


stations, it comes into notice that 
the microseismic waves originate 
definitely in the particularly lim- 
ited regions while the center of 
typhoon moves. The writer rec- 


ognized that the frequency of gen- 


eration is closely connected with a ; 
the slope of the continental shelf. _S wakiigin Equidepth line 


of 200m 
The equidepth lines of 200 meters 


drawn in thos 1g 5 Ss 
hose figures suggest 


i 


thermore clearly by Fig. 8. In Fig. 7 Frequency distribution of arrival dir- 


that fact, and it is indicated fur- 


ections of only microseismic waves having 


the typical Rayleigh form. (at Abuyama 
tion is supposed to be also associ- station ) 


this case the frequency of genera- 


ated with an angle of the arrival direction formed with the continental 


margin, besides a slope of the continental shelf, because the origin of 


80 
= 60 
be 
Sy) 
> 
& 40 

20 

0 

10 20 20 40 50 
D/ Sin @ 
(1) O- Yura station, @- Aso Laboratory 
Oo Japan Sea 
12 . ®@ Pacific Ocean 
10 


Frequency 
I a co 


nN 


i=) 


10 20 30 40 50 
D/Sin@ 


(2) Abuyama Observatory 


Fig. 8 Relationship between the frequency of microseismic waves and 
the ratio of D to Sin @. 
D : Distance from coasts to continental margins 
@ : Angle of arrival directions formed with continental margins 


microseisms is taken as the line source which is parallel with the continental 
margin. Therefore D/sin@ was adopted as the abscissa, in which D is the 
distance from the coast to the continenental margin and @ is the above 
mentioned angle. It may be readily understood from the above mentioned 
inference that there are minimum frequencies in the southerly direction 
in their distributions at the Aso station and also in the direction of the 
Wakayama district and the Noto Peninsula in the distributions at the 
Abuyama’ station. 

Let us consider why we observe frequently the microseismic waves 
coming from the regions possessing the continental shelves of steep slope. 


As the wave length of sea waves is comparable to the depth of the 


10 


continental shelf, the wave energy is dissipated as a result of friction by 
the oscillating motion of waves at the sea bottom, and in consequence 
the wave height reduces as the waves propagate shoreward. Accordingly 
as a decrease of the slope of shelf, the sea waves reduce in height in the 
vicinity of the coast, because of dissipating of the energy due to an increase 
of the distance of propagation, and then it results in the reduction of 
amplitudes of microseismic waves. The arrival direction may be thus not 


frequently observed by us in the direction of the gently sloping continental 
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Fig. 10 The amplitude of microseismic waves versus the arrival direction. 


shelf. If it is so, the amplitude of microseismic waves should be in pro- 
portion to the frequency of arrival directions. The relation between the 
amplitude and the arrival direction is shown in Fig. 9 for the microseismic 
waves coming from the Japan Sea and in Fig. 10 for those coming from 
the Pacific. The broken lines in those figures indicate the mean amplitudes. 
We may observe the upper limiting points plotted to be somewhat pro- 
portional to the frequency, but the fluctuation of the amplitude fad into 
insignificance in comparison with that of the frequency, especially for the 
mean amplitude. The writer could not thus recognize the close connection 
between the amplitude of microseismic waves and their frequency of gen- 
eration. From this fact it is hardly likely that microseismic waves with 
the energy enough to be propagated to the Abuyama station may be gen- 
erated at the regions close to the coastlines. The writer does not readily 
accept the consideration that the microseismic waves should be generated 
by the standing sea waves due to the interference of incident and reflected 
waves produced by a steep coast, because he wonders if the standing sea 
waves can so strongly generate elastic waves in the crust that the ground 
displacement is able to attain to a few microns after propagating the path 
of more than 100 km. 

According to the experimental study of Cooper and Longuet-Higgins 
(1951), the first order pressure variations by progressive water waves give 
considerable energy in the bed of tank, if the water depth is small as 
compared with half the wave-length. Therefore, the progressive sea waves 
must more strongly generate microseismic waves than the standing sea 


waves in sea beds close to coastlines, because the progressive waves, in 
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general, give more energy in phase in sea beds than the standing sea waves. 
The writer prefers the district at a distance from the coast, that is, 
rather near the continental margin, to the district close to the coastline as 
the origin of generation of microseisms. 
In conclusion steeper inclination of the sea bed must be accompanied 
by the condition to generate more frequently micorseismic waves, and at 
present we can but guess that the standing sea waves are frequently pro- 


duced under such a condition. 


6. Investigation of the wave period of microseisms 


coming from different arrival directions 


It is interesting and very significant to the investigation of microseisms 
to observe whether the microseismic waves with different periods come 
from different arrival directions. The observation to this purpose was 
carried out at the Abuyama station for the microseismic waves coming 
from the Japan Sea and the Pacific. 

When the progression of microseismic waves was previously investigated 
(Okano 1960), the ordinary seismograms in NS component were obtained 


together with the vector seismograms. Therefore the wave period could 
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Fig. 1] Relation of the period to the arrival direction of microseismic 
waves caused by seasonal winds in winter season. 
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be exactly read from the former and the arrival direction from the latter, 
and those relations are illustrated in Fig. 11. Although the period are 
considerably scattered centering around about 4.2 sec., on the whole there 
is little if any difference among those of waves coming from different 
directions. The writer supposed that those facts can be interpreted as 
follows. Since the seasonal wind stires widely up the Japan Sea in winter 
season, the disturbance source also covers a wide range of the Sea, and 
hence it is not so distant from the continental margin where the writer 
considers to be the origin of generation. Accordingly the wave period of 
microseisms is not so large because of the short distance of propagation 
of sea waves and it is not so different among the waves generated at the 


different districts. 


Aug. 20 1960 
20 pyhys™— 1929” 


ait 


Ee 
320 = 340 0 20 “40 60 60 100 120 140 160 160 200 220 240 


Direction ( dearee) 


Fig. 12 Relation of the period to the arrival direction of microseismic waves 
caused by typhoon No. 14. 


Fig. 12 shows the relation of the period to the arrival direction of 
microseismic waves when the typhoon no. 14 was situated as shown in 
Fig. 6. The difference of the periods among arrival directions is consid- 
erably noticed, excepting the waves propagated from the Japan Sea. Then 
we assume that the swells leave the end of the fetch shown in Fig. 6, 
propagate on the open sea and generate the microseismic waves in the 
vicinity of the continental margin. According to Sverdrup and Munk 
(1947), the relation of the wave period of swells at the end of the decay 


distance to that at the end of the fetch was indicated as follows. 


* of 1416n? Ar cee 


A: 27? 0'/p=6.35X 10-° 


14 


72 : Resistance coefficient applicable to wind 
(2.6 10-8 for wind velocities 5m/sec) 

: Density of air (1.25x 10~* g/em®) 

: Density of water (1.025 g/cm®) 


~ 


RF D> SD 


: Coefficient of energy partition (0.580) 
D : Distance of decay 


Based on the probable assumption that Ty is 6.0 sec., Tp was calcu- 
lated for each arrival direction, and the relation are indicated by the solid 
lines in Fig. 12, In this case the writer applied the Longuet-Higgins’s 
theory to the relation between the period of microseismic waves and that | 
of swells, according to which the former is half of the latter. Compara- 
tively good coincidence is recognized between the observation and the 
calculation. This is one of the evidences to show that the microseismic 
waves should be generated in the neighbourhood of the continental margin 
by the swells propagated from the disturbance source. The waves of the 
short period coming from the direction of 100-110 degrees may be the 3- 
second microseisms which are very frequently observed at Nagoya. (Na- 


tional Comittee for I.G.Y. 1959 and 1960). 
7. Summary 


l. The inference that microseismic waves are generated not in the 
neighbourhood of a low pressure center but at the regions near coasts, 
became quite obvious from the observations at the Aso and the Yura 
stations. 

2. The microseismic waves are more frequently generated at the dis- 
trict where the continental margin is more shoreward, namely, the slope 
of the continental shelf is steeper. 

3. The writer could not see the appreciable differences in mean 
amplitude among the waves coming from different directions. This fact 
suggests that microseismic waves should be generated at the districts not 
close to coastlines but rather near continental margins, because the swells 
are supposed to be alike in height everywhere in the vicinity of continental 
margins. 

4. The microseismic waves coming from different arrival directions 


are alike in mean period when the waves are generated by the seasonal 


LS 


wind in winter season. When the microseismic waves are generated by 
the typhoon, we may observe longer periods for the waves propagated from 


the district having longer distance of propagation of swells. 
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